arXiv:1509.02514v2 [astro-ph.SR] 29 Sep 2015 


Astronomy & Astrophysics manuscript no. Iindberg_arxiv2 ©ESO 2015 

September 30, 2015 


Probing the effects of external irradiation on low-mass protostars 

through unbiased line surveys^ 

J. E. Lindbergh’J. K. J0rgensen\ Y. Watanabe^, S. E. Bisschop\ N. Sakai^, and S. Yamamoto^ 

^ Centre for Star and Planet Formation, Niels Bohr Institute and Natural History Museum of Denmark, University of Copenhagen, 
0ster Voldgade 5-7, DK-1350 Copenhagen K, Denmark 

^ NASA Goddard Space Flight Center, Astrochemistry Laboratory, Mail Code 691, 8800 Greenbelt Road, Greenbelt, MD 20771, 
USA 

e-mail: johan.lindberg@nasa.gov 

^ Department of Physics, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan 


Received March 30, 2015; accepted September 7, 2015 


ABSTRACT 

Context. The envelopes of molecular gas around embedded low-mass protostars show different chemistries, which can be used to 
trace their formation history and physical conditions. The excitation conditions of some molecular species can also be used to trace 
these physical conditions, making it possible to constrain for instance sources of heating and excitation. 

Aims. We study the range of influence of an intermediate-mass Herbig Be protostar. We also study the effect of feedback from the 
environment on the chemical and physical properties of embedded protostars. 

Methods. We followed up on an earlier line survey of the Class O/I source R CrA IRS7B in the 0.8 mm window with an unbiased line 
survey of the same source in the 1.3 mm window using the Atacama Pathfinder Experiment (APEX) telescope. We also studied the 
excitation of the key species H 2 CO, CH 3 OH, and C-C 3 H 2 in a complete sample of the 18 embedded protostars in the Corona Australis 
star-forming region. Radiative transfer models were employed to establish abundances of the molecular species. 

Results. We detect line emission from 20 molecular species (32 including isotopologues) in the two surveys. The most complex 
species detected are CH 3 OH, CH 3 CCH, CH 3 CHO, and CH 3 CN (the latter two are only tentatively detected). CH 3 CN and several 
other complex organic molecules are significantly under-abundant in comparison with what is found towards hot corino protostars. 
The H 2 CO rotational temperatures of the sources in the region decrease with the distance to the Herbig Be star R CrA, whereas the 
C-C 3 H 2 temperatures remain constant across the star-forming region. 

Conclusions. The high H 2 CO temperatures observed towards objects close to R CrA suggest that this star has a sphere of influence 
of several 10 000 AU in which it increases the temperature of the molecular gas to 30-50 K through irradiation. The chemistry in the 
IRS7B envelope differs significantly from many other embedded proto stars, which could be an effect of the external irradiation from 
RCrA. 

Key words, stars: formation - ISM: individual objects (R CrA) - ISM: molecules - astrochemistry - radiative transfer 


1. Introduction 

The molecular composition of the envelopes of deeply embed¬ 
ded low-mass protostars can be used to trace their history and 
physical conditions. As a result of observations of these sources 
as well as of chemical modelling, the relevance and origin of 
many of these molecules are now relatively well understood 
(see e.g. Herbst & van Dishoeck 2009). The studied molecules 
include organic, silicon- and sulphur-bearing, and deuterated 
species (e.g. Blake et al. 1994; van Dishoeck et al. 1995; Schoier 
et al. 2002). Jprgensen et al. (2004) correlated the presence and 
abundance of nine molecular species, discovering strong corre¬ 
lations between some species and anti-correlations between oth¬ 
ers. This enabled constructing an empirical chemical network of 
molecular species. 


Based on observations with the Atacama Pathfinder Experiment 
(APEX) telescope. APEX is a collaboration between the Max Planck 
Institute for Radio Astronomy, the European Southern Observatory, and 
the Onsala Space Observatory. 

NASA Postdoctoral Program Fellow. 


Data from spectroscopy observations of molecular emission 
lines can also be used to trace the physical properties of proto- 
stellar envelopes. The rotational temperature of the molecular 
gas, which can be used as a proxy for the kinetic temperature 
of the gas given certain assumptions, can be measured by exci¬ 
tation diagrams of the observed molecules (see e.g. Goldsmith 
& Danger 1999; Jprgensen et al. 2005b). By comparing the rota¬ 
tional temperatures of different molecular species, insights about 
their distribution and formation paths can be gained. Isotopical 
ratios can also be used as tracers of formation history. 

During the past decade, the existence of an inner (R < 
100 AU) hot region in the envelopes of deeply embedded low- 
mass protostars with a chemistry similar to that of high-mass 
hot cores has been proposed. These dense (n > 10^ cm“^) re¬ 
gions are referred to as hot corinos (Ceccarelli 2004), and they 
are characterised by the presence and relatively high abundance 
of complex organic molecules such as CH 3 OH, CH 3 OCH 3 , 
and CH 3 OCHO (Bottinelli et al. 2004a, 2007; Jprgensen et al. 
2005a). High angular resolution interferometric observations 
show that complex organics are often concentrated in the in¬ 
ner < 1" (R < 50-100 AU) of low-mass protostars (e.g. Bot- 
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tinelli et al. 2004b; Kuan et al. 2004; J 0 rgensen et al. 2005a; 
Bisschop et al. 2008; Maury et al. 2014) - most likely reflecting 
that they evaporate from the icy dust grains at high temperatures 
{T > 100 K). IRAS 16293-2422 is considered as the prototyp¬ 
ical hot corino, but at least three more such sources have been 
discovered. 

Sakai et al. (2008, 2009b) found that the chemistry towards 
two deeply embedded low-mass protostellar sources was dis¬ 
tinctly different from that of the hot corinos. The two sources 
L1527 and IRAS 15398-3359 are dubbed warm carbon-chain 
chemistry (WCCC) sources. They are characterised by compa¬ 
rably high abundances of long carbon-chain molecules such as 
HC 3 N, C 4 H, and C 6 H, and at the same time by low or mod¬ 
erate abundances of complex organic species. The difference 
between the chemical properties observed towards hot corinos 
and WCCC sources could be an effect of different collapse 
timescales (Sakai et al. 2009b; Sakai & Yamamoto 2013). This 
mechanism has been verified by chemical modelling (Hassel 
et al. 2008). 

The presence of complex organic molecules in low-mass star 
formation is, however, not unique to the hot inner regions of 
deeply embedded protostars. The spectral signatures of several 
complex organic molecules have recently been detected towards 
the prestellar core L1689B, with temperatures ~ 10 K (Bacmann 
et al. 2012). Complex organic molecules have also been found 
towards protostellar outflows such as that of LI 157, with simi¬ 
larly low temperatures (Arce et al. 2008). 

The transitional Class O/I protostar R CrA IRS7B is located 
in the R CrA cloud (NGC 6729) within the Corona Australis 
(CrA) star-forming region. At a distance of 130 pc (Neuhauser & 
Forbrich 2008), this is one of the nearest low-mass star-forming 
regions, making it ideal for the study of the envelope and core 
chemistry of the individual sources by single-dish (sub)mm ob¬ 
servations. The population of young stellar objects (YSOs) in 
the region was inventoried through multi-wavelength studies by 
Peterson et al. (2011), who identified 116 YSOs and classified 
them according to their spectral slope following the definitions 
of Greene et al. (1994). Out of these sources, 14 were found to be 
Class I sources or younger, and 5 were found to be flat-spectrum 
sources (transitional sources between Class I and Class II). Most 
of these deeply embedded sources are concentrated in the R CrA 
cloud, within a few 1000 AU from the Herbig Be star R CrA. 
Lindberg & Jprgensen (2012) found large-scale H 2 CO emission 
in the region, exhibiting gas temperatures > 30 K. The high tem¬ 
peratures, which cannot be explained by internal radiation from 
the low-mass stars, were suggested to be a result of external ir¬ 
radiation from R CrA. 

Watanabe et al. (2012) presented an unbiased line survey of 
R CrA IRS7B in the 0.8 mm window and a few scans in the 
0.7 mm window, using the Atacama Submillimeter Telescope 
Experiment (ASTE) 10 m telescope. These observations showed 
enhancements of the radicals CN and C 2 H, which trace photon- 
dominated regions (PDRs). Eurthermore, neither any complex 
organic molecules characteristic for hot corino chemistry (except 
the ubiquitous CH 3 OH) nor any long carbon chains characteris¬ 
tic for WCCC were detected. High-resolution Atacama Large 
Millimeter/submillimeter Array (ALMA) observations towards 
R CrA IRS7B also indicate an absence or unusually low abun¬ 
dances of the complex organic species CH 3 CN and CH 3 OCH 3 
(Lindberg et al. 2014a). 

This paper complements the work by Watanabe et al. (2012) 
by adding an unbiased line survey in the 1.3 mm window, which 
increases the number of detected molecular species as well as 
the number of molecular transitions in many of the excitation 


diagrams, and thus also the range of in those diagrams. We 
also investigate other similar protostellar sources in Corona Aus¬ 
tralis to study the spatial variations of the chemical and physical 
properties in the star-forming region. This paper finally uses ra¬ 
diative transfer modelling of the emission from several of the 
species to establish molecular abundances for comparison with 
other sources. 

Section 2 describes the observational methods and data re¬ 
duction strategy. Section 3 describes the spectra resulting from 
the observations, and Sect. 4 discusses rotational diagrams, ra¬ 
diative transfer models, and isotopologue ratios of many of the 
observed molecules. Section 5 discusses the large-scale heating 
in the CrA region and the chemistry of the envelopes. Section 6 
lists our conclusions. 


2, Observations 

2.1. APEX unbiased line survey of R CrA IRS7B 

The observations were carried out with the Atacama Pathfinder 
Experiment 12 m telescope (APEX; Giisten et al. 2006) in 
position-switching mode in August and October 2010. The fre¬ 
quency range between 217.9 GHz and 245.5 GHz was covered 
with 17 spectral setups on the SHeEI (Swedish Heterodyne Ea- 
cility Instrument) receiver APEX-1 (Vassilev et al. 2008). This 
survey was unbiased in the sense that it covered a wide spectral 
range without targeting any particular molecular species. In ad¬ 
dition, we used parts of the APEX-1 data that were observed by 
Lindberg & Jprgensen (2012) in April 2010 to cover the range 
between 217.2 GHz and 217.9 GHz. All observations were cen¬ 
tred at the sub-mm/IR point source (Class O/I YSO) IRS7B in 
the R CrA cloud (see Table 1 for coordinates). The two sets of 
observations have an overlap around 218 GHz, which we used to 
investigate the calibration accuracy. The small differences found 
in the line intensities can all be explained with a calibration un¬ 
certainty < 10 %. 

Eor all new scans reported in this work, rms levels between 
21 and 35 mK channel"^ were reached with velocity resolutions 
ranging between 0.17 km s“^ and 0.15 km s“^ (122 kHz), cor¬ 
responding to rms levels of 8-14 mK (km s“^)“^ The scan 
from Lindberg & Jprgensen (2012), which here was used to 
cover 217.2 GHz-217.9 GHz, is not as deep as the new data 
presented here. It has an rms level of 84 mK channel" \ or 
35 mK (km s"^)"^ 

Because of unstable quasi-sinusoidal baselines in the ob¬ 
served spectra (see Vassilev et al. 2008), we had to perform rela¬ 
tively advanced baseline fitting after carefully identifying chan¬ 
nels containing any spectral line emission. The exact method 
used for this baseline fitting is described in Appendix A. The cal¬ 
culation and subtraction of baselines were performed using our 
own scripts as described in the appendix, but Gaussian line fit¬ 
ting and calculations of line strengths were afterwards performed 
using the GILDAS CLASS packaged To account for some 
smaller features left after the baseline subtractions, a zeroth- 
order baseline was subtracted locally around all fitted lines be¬ 
fore the extraction of line parameters. 

All line intensities are given in the Ymb (main beam bright¬ 
ness temperature) scale, which means that the antenna temper¬ 
atures are corrected with the APEX main beam efficiency 
= 0.75. The beam sizes of the observations range between 

^ GILDAS CLASS (Continuum and Line Analysis Single-dish Soft¬ 
ware) is developed by the IRAM institute, Grenoble, France, and can 
be downloaded from http: //www. iram. fr/IRAMFR/GILDAS. 
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Table 1. Observed sources and the rms of the observations in the APEX surveys of CrA. See Table 2 for spectral window parameters. 


Source name^ 

RA 

(J2000.0) 

Dec 

(J2000.0) 

YSO Class*’ 

rms [mK (km s 
1.4 mm 0.9 mm 

i)-i] 

0.8 mm 

IRS7B 

19 01 56.40 

-36 57 28.1 

I 

IE 

13 

40 

CrA-46 

18 55 56.32 

-37 00 07.1 

Elat 

19 



CrA-3 

18 59 43.92 

-37 04 01.1 

Elat 

19 

32 


CrA-5 

19 00 15.55 

-36 57 57.7 

I 

10 

21 

27 

LS-RCrAl 

19 01 33.56 

-37 00 30.3 

Elat 

20 

22 


Haas4 

19 01 40.67 

-36 56 04.9 

Elat 

16 

20 


IRS2 

19 0141.56 

-36 58 31.2 

I 

7 

14 


IRS5A 

19 01 48.03 

-36 57 22.2 

I 

11 

18 

36 

IRS5N 

19 01 48.46 

-36 57 14.7 

I 

11 

15 

38 

IRSl 

19 01 50.68 

-36 58 09.7 

I 

12 

20 

34 

IRS7A 

19 01 55.32 

-36 57 21.9 

I 


14 

38 

CrA-24 

19 01 55.60 

-36 56 51.1 

I 

20 

20 

36 

SMM 2 

19 01 58.54 

-36 57 08.5 

I 

20 

15 

36 

CX042 

19 02 01.96 

-36 54 00.0 

I 

8 

17 


CrA-44^ 

19 02 58.67 

-37 07 35.9 

I 

12 

19 


CrA-33 

19 03 01.03 

-37 07 53.4 

Elat 

19 

16 


VV CrA 

19 03 06.80 

-37 12 49.1 

I 

9 

12 


CrA-37 

19 03 55.24 

-37 09 35.9 

I 

20 

15 



Notes. We use the notation of Peterson et al. (2011) for all sources except IRS5A, IRS5N, and SMM 2, for which we use the notation of 
Lindberg & Jprgensen (2012) and Lindberg et al. (2014b). These are referred to as CrA-19, CrA-20, and CrA-43, respectively, by Peterson et al. 
(2011). Assigned from the spectral index a (Peterson et al. 2011). Class I sources are Class I or younger. Only partially covered by the 
IRS7B line survey. The observations do not cover 216.1-217.2 GHz, and the 217.2 GHz-217.9 GHz portion is covered only at a lower S/N level 
(37 mK (km Also known as IRAS 32. 


26" and 29" depending on the observed frequency. All errors 
given in this work are at the Icr level if not otherwise stated. 

As a complement to the APEX observations, we use the 
spectral lines detected in the ASTE 0.8 mm and 0.7 mm ob¬ 
servations of Watanabe et al. (2012). This enables a more accu¬ 
rate excitation analysis of the various molecular species thanks 
to the larger range of E’u in the used transitions. The beam size 
of these ASTE observations ranges between 21" and 23" for the 
0.8 mm window, and between 16" and 17" for the few lines 
in the 0.7 mm window. Eor the 0.8 mm window in particular, 
the beam size is quite similar to that of the APEX observations, 
meaning that the two sets of observations cover roughly the same 
spatial region. Erom previous SMA observations we know that 
the spatial extent of the molecular line emission is similar to or 
larger than these beams (Lindberg & Jprgensen 2012). Eor this 
reason, and to be consistent with Watanabe et al. (2012), we thus 
calculated beam-averaged column densities. When combining 
data points with different beam sizes in excitation diagrams, this 
means that we assumed that the beam-averaged column densities 
are similar in both beams, or in other words, that the emission is 
extended and uniform. 


2.2. Unbiased survey of H2CO, CH3OH, and C-C3H2 in CrA 

In addition to the unbiased line survey of IRS7B discussed 
above, IRS7B and 17 other protostellar sources within the 
Corona Australis star-forming region were observed in more nar¬ 
row frequency bands using the APEX telescope. The selected 
sample consists of all sources defined as Class I sources (or 
younger) or flat-spectrum sources in the survey of YSOs in CrA 
conducted by Peterson et al. (2011). Only the faint Class I source 
CX034 was excluded because of its close proximity to the two 


Table 2. Observational parameters of the APEX observations (this 
work; IRS7B line survey and CrA source survey) and the ASTE line 
survey of IRS7B (Watanabe et al. 2012). 


Window 

Erequency range 

Beam size 


[GHz] 

["] 

[AU] 

IRS7B survey 

217.2-245.5^ 

26-29 

3300-3800 

1.4 mm 

216.1-220.1 

29 

3700 

0.9 mm 

336.9-340.9 

19 

2400 

0.8 mm 

361.8-365.8*’ 

17 

2200 

ASTE 0.8 mm 

332.0-364.0 

21-23 

2700-3000 

ASTE 0.7 mm 

435.0-437.0 

17 

2300 

ASTE 0.7 mm 

459.2-461.2 

16 

2100 


Notes. The 217.2 GHz-217.9 GHz portion is only covered at a low 
S/N level. In most sources only 361.8-364.3 GHz, due to high noise 
levels and baseline issues in the upper sideband. 


much stronger sources IRS7A and IRS7B, which would cause 
confusion. The observed protostellar sources are shown on a 
map in Eig. 1, and their positional data are listed in Table 1. The 
observations were performed in up to three frequency setups: 
216.1-220.1 GHz (1.4 mm), 336.9-340.9 GHz (0.9 mm), and 
361.8-365.8 GHz (0.8 mm), see also Table 2. The frequency se¬ 
tups were chosen to cover a wide range of H 2 CO, CH 3 OH, 
and C-C 3 H 2 lines. The 1.4 mm window was observed with the 
SHeEI APEX-1 receiver and the other two windows with the 
SHeEI APEX-2 receiver (Vassilev et al. 2008). Observations in 
the higher-frequency windows were only executed for sources 
with a significant number of detected lines (see Table 1 for de- 
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Fig. 1. Overview of the observed sources in the H 2 CO, CH 3 OH, and C-C 3 H 2 survey. The circles show the APEX beam size in the 1.4 mm 
window. Red circles denote Class O/I objects and blue circles are flat-spectrum sources. CrA-46 is located ~ 40' (~ 1.5 pc) west of the right edge 
of the plot. The inset shows a zoom-in of the region around R CrA enlarged by a factor ~ 2. The greyscale background is a Spitzer 4.5 [im image. 
(A colour version of this plot is available in the online journal.) 


tails on observed spectral windows). The observations were car¬ 
ried out in September and November 2011. This survey is unbi¬ 
ased in the sense that it covers (nearly) all embedded sources in 
CrA, albeit only in a limited number of spectral setups, targeting 
certain molecular species. 

All line intensities are given on the scale. The APEX 
main beam efficiency = 0-75 for the 1.4 mm window and 
= 0.73 for the 0.9 mm and 0.8 mm windows. The un¬ 
smoothed spectra have 76 kHz channel widths (corresponding 
to 0.10 km s“^ for the 1.4 mm window, 0.067 km s“^ for the 
0.9 mm window, and 0.063 km s“^ for the 0.8 mm window). 
The rms levels when smoothed to 1 km s"^ channels are listed 
in Table 1. In the 1.4 mm window, the rms level is between 
7 mK (km s“^)“^ and 20 mK (km s“^)“\ which can be compared 
to a noise level of 11 mK (km s“^)“^ for the 1.4 mm window in 
the IRS7B survey (above). 

Like in the observations described in Sect. 2.1, the spectra in 
this survey also have unstable baselines, and the same baseline 
corrections as were used for the IRS7B data were applied (see 
Appendix A). The 0.8 mm band was particularly affected by 
unstable baselines with both wide and narrow features. In the 
lower sideband of the 0.8 mm observations these effects could 
mostly be accounted for, but in the upper sideband the features 
were too strong to make proper use of the spectra except for 
in smaller stable chunks. Thus, only a few spectral lines were 
identified in this sideband, and then only in the sources where 
the strongest line emission was found. 

3. Results 

3.1. IRS7B line survey 

We detect 102 spectral lines in the IRS7B line survey between 
217.2 GHz and 245.5 GHz. We successfully identify 87 of the 
detected lines, and the identifications correspond to 16 molecu¬ 
lar species (25 if counting isotopologues separately). The online 


database Splatalogue^ was used for line identification. The de¬ 
tected spectral lines and their parameters are listed in Table B.l 
in Appendix B. A compressed version of the APEX spectrum of 
the IRS7B observations can be found in Eig. 2, and a full spec¬ 
trum is available in Appendix C (Eigs. C.1-C.8). 

In Eig. 3, we compare the line emission observed towards 
IRS7B in a window around 227 GHz with the emission detected 
towards a typical high-mass hot core and a typical low-mass hot 
corino source. This simplistic comparison shows that the hot 
corino bears more resemblance to the hot core than to IRS7B. 

Many of the species detected towards IRS7B were previously 
detected in the ASTE survey (Watanabe et al. 2012), but we de¬ 
tect six molecular species not detected in that survey: HNCO, 
H 2 CCO, CH 3 CCH, CH 3 CHO, CH 3 CN, and HC 3 N. HC 3 N was 
previously detected in SMA/APEX observations (Lindberg & 
Jprgensen 2012), but the five other species have, to the best of 
our knowledge, not previously been detected towards IRS7B. 
The CH 3 CN and CH 3 CHO detections are only tentative: Eor 
CH 3 CN, only two faint spectral lines are seen, and the two lines 
at 239.133 GHz and 239.138 GHz, expected to be of similar in¬ 
tensity, are not detected. Eor CH 3 CHO, three lines are detected 
at 3-6cr levels, but three lines at 219.820 GHz, 232.691 GHz, 
and 235.684 GHz, expected to be at least as strong, are non- 
detected. We also detect some isotopologues not seen in the 
ASTE survey: ^^CN, C^^N, ^^CO, C^^O, ^^CS, C^^S, and DNC. 

We find that all identified spectral lines have LSR velocities 
in the range 5.0-6 .6 km s“^ (median 5.7 km s“^). The line widths 
show a somewhat larger variation, with the lowest values around 
1 km s“^ for species like ^^CN, C-C 3 H 2 and CH 3 CCH, and the 
highest values above 6 km s“^ for CS and CO. The median line 
width for the non-blended lines to which a Gaussian could be 
fitted is 2.0 km s“^ 

Of the 102 detected lines, 15 are unidentified. A handful of 
these are on the 3-5cr level, and could thus be flukes, but a few 

^ http://splatalogue.net/ 
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Rest Frequency [GHz] 


Fig. 2. 1.3 mm spectrum of IRS7B, cut at 10.5 K (top panel) and 0.5 K (bottom panel). A high-resolution spectrum can be found in Appendix C. 


U lines are strong. In Appendix E we demonstrate our attempts 
to identify these U lines. 

3 . 2 . H2CO, CH3OH, and C-C3H2 survey in CrA 

Of the 18 targeted sources in this survey, only CrA-46 shows no 
spectral line emission at all. Towards an additional four sources, 
only CO isotopologues were detected, but towards the remain¬ 
ing 13 sources several different molecular species were detected. 
The species that were detected towards at least one of the sources 
in this source survey are H 2 CO, H^^CO, CH 3 OH, C-C 3 H 2 , C^^O, 
C^^O, DCO+, HCi^O+, SiO, SO, ^^SO, CN, ^^CN, C^^N, 

C 2 D, DCN, HNC, and H 2 CS. The detected spectral lines and 
their parameters are listed in Tables B.2-B.22 in Appendix B, 
except for the IRS7B parameters, which are given together with 
the results of the IRS7B line survey in Table B.l. The spectra 
can be found in Figs. D.1-D.41 in Appendix D. 

With the exception of SiO, which traces shocked gas at high 
absolute velocities relative to the protostar, all spectral lines have 
LSR velocities in the range 4. 6 - 6.6 km s“^ towards all sources 
where spectral lines are detected (CrA-24 and SMM 2 have two 
velocity components, however, which is discussed below). This 
confirms the conclusion of Peterson et al. (2011) that all sources 
where we detect line emission are situated in the same star¬ 
forming region, none of them representing background sources. 

In CrA-24 and SMM 2, the H 2 CO, CH 3 OH, SO, and SiO 
spectral lines show two velocity components: one low-velocity 
component with LSR velocities in the range 1.2-3.1 km s“^ (out¬ 
flow component) and one component in agreement with the other 
sources, with LSR velocities in the range 4.6-6.4 km s“^ (on- 
source component). The SiO emission is strongest in the low- 
velocity component, while the other species are stronger in the 


on-source component. The low-velocity component is probably 
an effect of a spatial overlap with a large-scale (~ 10000 AU) 
outflow originating in IRS7B, extending in a north-eastward di¬ 
rection. This outflow has been detected in SEST mapping obser¬ 
vations of the CH 3 OH 2-2 ^ 1 - 1 , E and 2o ^ lo, A-r transi¬ 
tions at 96.7 GHz (Miettinen et al., in prep.). Traces of this out¬ 
flow can also be seen in ALMA observations of CH 3 OH lines at 
340 GHz (Lindberg et al. 2014a). For the CrA-24 and SMM 2 
data, we separated the line emission of the two velocity compo¬ 
nents by fitting two Gaussians to each spectral line. The total line 
intensities, the on-source components, and the outflow compo¬ 
nents are reported in three separate tables for each of these two 
sources in Appendix B. The lines of all species except H 2 CO, 
CH 3 OH, SO, and SiO show only one velocity component, which 
is consistent with the on-source component. These lines are only 
listed in the tables with total line intensities. 

We note that most sources with a large number of line de¬ 
tections are found in the R CrA cloud. CrA-44 and VV CrA 
also show large numbers of line detections, while the remaining 
sources far from R CrA show relatively few line detections. 


4. Analysis 

4.1. Rotational diagram analysis - IRS7B 

We performed a rotational diagram analysis (see e.g. Goldsmith 
& Langer 1999) of several of the detected molecular species. 
This analysis requires the assumptions that the gas is in local 
thermal equilibrium (LTE) and that all lines are optically thin. In 
addition to the APEX data, we included ASTE data previously 
reported by Watanabe et al. (2012) in the fits. We performed this 
analysis for all the molecular species detected towards IRS7B for 
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Fig. 3. Single-dish spectra around 227 GHz of the hot core G327.3-0.6 
(top; observed with APEX; Bisschop et al. 2013; Bisschop et ah, in 
prep.), the hot corino IRAS 16293-2422 (middle; observed with IRAM 
30 m; Caux et al. 2011), and R CrA IRS7B (bottom; this work). The 
labels in the upper plot correspond well to the lines also seen in the 
middle plot. The colours denote oxygen-bearing (purple), nitrogen¬ 
bearing (blue), sulphur-bearing (yellow), and typcial WCCC and PDR 
(red) species. (A colour version of this plot is available in the online 
journal.) 


which at least two transitions with significantly different were 
detected (a minimum difference of 15 K was chosen to avoid 
making rotational diagram fits with large errors). This includes 
H 2 CO, D 2 CO, H^^CO, CH 3 OH, CH 3 CCH, H 2 CS, SO, SO 2 , CS, 
C^^S, CN, DCN, HC 3 N, C 2 H, C-C 3 H 2 , and CH 3 CHO. To estab- 
fish approximate column densities, we also performed fits for 
H 2 CCO, HNCO, and CH 3 CN assuming a fixed temperature of 
30 + 5 K. We also fixed the rotational temperature to 30 + 5 K 
for CH 3 CHO because the three detected lines are relatively faint. 
The temperature 30 K was chosen to match the rotational tem¬ 
peratures measured for other organic species. We excluded CO 
and its isotopologues because their lines are probably optically 
thick. We used the CDMS (Muller et al. 2001) and JPL (Pickett 
et al. 1998) databases to retrieve molecular spectroscopy data. 

From previous interferometric observations we know that the 
line emission in the region is very extended in comparison with 


the APEX and ASTE beams (Lindberg & Jprgensen 2012), and 
we therefore calculated beam-averaged column densities. This 
requires the assumption that the column densities are relatively 
uniform across the range of beam sizes. Since the 0.9 mm and 
0.8 mm band APEX observations of IRS7B partially coincide 
with the ASTE survey, we tested this assumption by comparing 
line strengths of lines covered in both surveys. We find that lines 
of H 2 CO, H 2 CS, HNC, DCN, and most CN lines are within er¬ 
rors between the two measurement sets, but lines of C^^O, C^'^S, 
and CH 3 OH are 1.5-2.3 times stronger in the ASTE data. This 
might be explained by the different spatial distribution of these 
molecules in combination with the difference in the APEX and 
ASTE beam sizes - the interferometric observations of Lind¬ 
berg & Jprgensen (2012) show that the CH 3 OH emission peaks 
~ 15" south of IRS7B, and this emission probably contributes 
more to the signal in the larger ASTE beam than in the APEX 
beam. When lines were available from both datasets, the ASTE 
data were used instead of the APEX data because the S/N ra¬ 
tio is higher in the ASTE data and because the APEX 1.3 mm 
beam size is more similar to the ASTE 0.8-0.9 mm beam sizes 
than the APEX 0.8-0.9 mm beam sizes. Three H 2 CO lines de¬ 
tected in the APEX 0.8 mm observations with frequencies just 
above the upper boundary of the ASTE survey (364 GHz) were 
included, however. To account for the calibration uncertainty of 
the observations, we introduced a 10 % calibration error in addi¬ 
tion to the rms noise. 

Eor the H 2 CO and H^^CO lines we assumed an ortho-to-para 
ratio of 1.6 (see Dickens & Irvine 1999; Jprgensen et al. 2005b), 
for D 2 CO we assumed the statistical value of 2 , and for C-C 3 H 2 , 
H 2 CS, and H 2 CCO we assumed ortho-to-para ratios of 3 (Lucas 
& Liszt 2000; Minh et al. 1991; Ohishi et al. 1991). Eor the 
molecules with A and E states (CH 3 OH, CH 3 CCH, CH 3 CN, and 
CH 3 CHO), we used the statistical A/E ratios as provided by the 
CDMS database. 

Mangum & Wootten (1993) showed that line ratios of H 2 CO 
transitions with the same /u-level are excellent probes of the ki¬ 
netic temperature of the molecular gas. However, line ratios of 
transitions involving different /u-levels also strongly depend on 
the density ^(H 2 ). This is further discussed in Appendix E. To 
account for this effect, we calculated rotational temperatures for 
each /u-level separately. We then used the weighted average of 
these two as the H 2 CO temperature, and found the column den¬ 
sity A^(H 2 C 0 ) and the molecular density ^(H 2 ) by comparing the 
measured 3 o 3 ^ 2 o 2 and 5 o 5 ^ 4 o 4 line strengths and RADEX 
estimates of these transitions (see Appendix E for a thorough de¬ 
scription of the method used). The H 2 number density towards 
IRS7B is found to be (8.9 + 1.0) x 10^ cm“^. 

The same problem probably applies to H^^CO and D 2 CO, 
but the low number of detected lines, their low S/N levels, and 
the lack of collisional data for these species make it difficult to 
perform fits on the separate /u-level transitions. The fits to these 
two isotopologues are thus highly uncertain. 

The results of all the IRS7B rotational diagram fits are listed 
in Table 3, and the rotational diagrams are shown in Eig. 4. The 
rotational temperatures of the different molecules are found to 
vary between 10 K and 40 K. Erom their rotational temperatures, 
the molecules can be divided into two groups: on one hand the 
nitrogen-bearing species (i.e. cyanides) and unsaturated hydro¬ 
carbon molecules at 10-17 K, and on the other hand, the oxygen¬ 
bearing organic molecules and sulphur-bearing species at 19- 
40 K. H^^CO, however, has a lower rotational temperature, but 
the fit is made with only three lines, involving transitions with 
different J^. 
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Fig. 4 . Rotational diagrams of molecular species detected towards IRS7B. Red square data points are 1.3 mm lines observed with APEX, 
blue circle data points are 0.8 mm lines observed with APEX (this work); green diamond data points are 0.8 mm lines observed with ASTE, and 
magenta hexagon data points are 0.7 mm lines observed with ASTE (Watanabe et al. 2012). Upper limit data points are semi-transparent. For 
H 2 CO, H 2 ^C 0 , D 2 CO, C-C 3 H 2 , and H 2 CS, the ortho lines are filled in the left half, the para lines are filled in the right half, and lines that are blends 
of ortho transitions and para transitions are completely filled. For CH 3 OH and CH 3 CCH, the A lines are filled in the top half and the E lines are 
filled in the bottom half. (A colour version of this plot is available in the online journal.) 


The rotational diagram of CH 3 OH shows a considerable 
scatter (see Fig. 4), which like in the case of H 2 CO suggests non- 
LTE conditions. To be able to accurately evaluate the tempera¬ 


ture and column density of this species, we employed non-LTE 
RADEX models. Erom interferometric observations, we know 
that H 2 CO and CH 3 OH are present in the same spatial regions 
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Table 3 . Rotational diagram parameters for IRS7B. 


Molecule 

Ttot 

[K] 

A^rot 

[cm- 2 ] 

HzCO'* 

40.3 + 2.9 

(1.0 ± 0.1) X 10'^ 

D2CO 

31.2 + 4.2 

(1.8±0.3)x 10'2 

HfCO 

13.6+1.8 

(3.4± 1.3)x 10‘2 

CH3OH 

27.6 ± 0.9 

(1.6±0.1)x 10'“^ 

CH3CCH 

32.9 + 2.1 

(1.1 ±0.2)x 10'“^ 

H2CS 

23.2 + 2.1 

(8.4±2.1)x 10'2 

SO 

28.9 + 1.9 

(7.1 ± l.l)x 10‘3 

S 02 

38.0 + 2.5 

(1.8±0.2)x 10‘3 

cs 

23.8 ± 2.6 

(7.3± 1.8)x 10'3 


18.5 ± 1.9 

(6.4± 1.7)x 10'^ 

CN 

9.9 + 0.3 

(1.8±0.1)x 10‘4 

DCN 

10.9 + 0.9 

(2.6 ± 0.6) X 10‘2 

HC3N 

15.6 ± 1.9 

(6.1 ±6.6) X 10*4 

C2H 

16.9 ±2.8 

(2.8± 1.2)x 10'“^ 

C-C3H2 

16.9 ± 0.7 

(8.5 ± 0.9) X 10'2 

HNCO 

30*’ 

(5.1 ± 1.9)x 10‘2 

CH3CN 

30'’ 

(8.7 ± 3.9) X 10" 

CH3CHO 

30'’ 

(1.6±0.7)x 10'^ 


Notes. The H 2 CO fit is a weighted average of the /u = 3 and /u = 5 
fits. The two fits were also used to estimate the H 2 number density, 
w(H 2 ) = (8.9 + 1.0) X 10^ cm"^. For these species, the Eyy of the ob¬ 
served lines are in too narrow a range for the rotational temperature to 
be accurately calculated, or all detected lines are of a very low S/N level. 
We assumed a temperature of 30 + 5 K (based on the rotational temper¬ 
atures of similar species) to be able to estimate the column densities of 
these molecules. 

around IRS7B (Lindberg & J 0 rgensen 2012), and we therefore 
assumed the H 2 density calculated from the H 2 CO excitation 
analysis (^(H 2 ) = 8.9 x 10^ cm ^). 

Our CH 3 OH models show that with the expected physical 
properties, a scatter in the CH 3 OH rotational diagram similar to 
what we observe should appear due to non-LTE excitation of 
the molecule. We performed a least-;^^ fit with the temperature 
and column density as free parameters. The errors used in the 

calculations are a combination of the rms error and a 10 % 
calibration uncertainty. We find that the lowest is achieved 
for r = 28 K and A^(CH 30 H) = 2.1 x lO^^ ^^^-2 best- 
fit temperature and column density are similar to the rotational 
diagram solution (see below). Figure 5 shows a contour plot of 
the reduced-;^^ values of fit. 

The median line width of the nitrogen-bearing and hydrocar¬ 
bon species ( 1.6 km s“^) is slightly lower than the median line 
width of the organic and sulphur-bearing species ( 2.1 km s“^). 
The median LSR velocity is not significantly different between 
the two groups of molecules. The differences in temperature 
and velocity dispersion suggest that the two groups of molecules 
have different spatial origins. Interestingly, molecules from the 
two groups were identified as having two different spatial dis¬ 
tributions in interferometric spectral line mapping of the region 
(Lindberg & Jprgensen 2012): the warmer molecules H 2 CO 
and CH 3 OH were found in two ridges north and south of the 
YSOs in the R CrA cloud, whereas the cooler molecules DCN, 
CN, C-C 3 H 2 , and HC 3 N were found to peak in an east-west 
band around the protostars and/or in a faint outflow extending 
northwards from IRS7B. CH 3 CCH does not exhibit the expected 
properties - despite being an unsaturated hydrocarbon and hav- 



Fig. 5 . Reduced-;^^^ fit to the CH 3 OH lines observed towards IRS7B 
assuming w(H 2 ) = 8.9 X 10^ cm"^. The contours indicate Icr, 2 cr, and 
3 cr certainty. The star symbol shows the best fit. 

ing a small velocity dispersion, it does have a rotational temper¬ 
ature similar to the second group of molecules, so its nature in 
relation to the other species remains unclear. 

4.2. Rotational diagram analysis - other sources 

Rotational diagrams of H 2 CO, CH 3 OH, and C-C 3 H 2 were pro¬ 
duced for the sources included in the CrA survey where at least 
two lines of the respective molecule were detected. For sources 
with only one detected line of a given species, upper limits of 
the rotational temperature and lower limits of the column den¬ 
sity were calculated. As above, we assumed ortho-to-para ra¬ 
tios of 1.6 for H 2 CO and 3 for C-C 3 H 2 , and an A/E ratio of 1 
for CH 3 OH. For the IRS7B rotational diagrams in this section 
we used only the APEX observations, and only those parts of 
the IRS7B survey that are also covered by the source survey, 
to avoid any bias between IRS7B and the other sources, even 
though more lines are available for IRS7B in the full line survey 
and the ASTE data. Thus, the rotational diagram parameters for 
IRS7B reported in this section are not identical to those in the 
previous section. We used the same calibration uncertainty of 
10% as in the rotational diagrams for IRS7B. 

Since no mm line interferometry data are available for most 
of the sources in the sample, we cannot determine the spatial ex¬ 
tent of the molecular line emission. We therefore used the same 
approach as for IRS7B, calculating the beam-averaged column 
densities. 

As discussed in Sect. 4.1, H 2 CO rotational diagrams must 
treat transitions involving different /u-levels separately. In the 
sources where lines from both the = 3 and = 5 levels 
are detected, we thus calculated two rotational temperatures and 
used the weighted average of these. We then estimated the col¬ 
umn density and the number density with the method described 
in Appendix F. We find that towards the sources where we could 
measure the H 2 density, it is between 4 x 10^ cm“^ and 10^ cm“^, 
except for towards the outflow component of SMM 2, where the 
density is lower. It is difficult to constrain the densities and col¬ 
umn densities for the outflow components of SMM 2 and CrA- 
24, since the densities are low and the column densities are high, 
which enhances both the optical depth effects and the non-LTE 
effects. 
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Fig. 6. H 2 CO rotational diagrams of the sources in the CrA survey where at least two H 2 CO lines were detected. The ortho lines are filled in 
the left half and the para lines are filled in the right half. Red square data points are 1.3 mm lines and blue circle data points are 0.8 mm lines, all 
observed with APEX. (A colour version of this plot is available in the online journal.) 


The rotational diagrams are shown in Figs. 6-8 and the fitted 
parameters are listed in Tables 4-5. In Fig. 9, the derived H 2 CO, 
CH 3 OH, and C-C 3 H 2 rotational temperatures are plotted versus 
the estimated distance to R CrA (the statistically most likely 
distance assuming a uniform three-dimensional spherical cloud 
with appropriate error bars are used to illustrate the uncertainty 
of projected distances). We find that the H 2 CO temperatures 
show a large spread throughout the cloud (19-74 K), the CH 3 OH 
temperatures show a smaller spread (12-46 K), while the C-C 3 H 2 
temperatures are uniform around 9-16 K. In the sources CrA- 
24 and SMM 2, the two velocity components of the H 2 CO and 
CH 3 OH spectral lines (corresponding to on-source and outflow 
emission, see Sect. 3.2) are treated separately. The H 2 CO tem¬ 
peratures of the outflow components are signiflcantly higher than 
the H 2 CO temperatures of the on-source components. This is 
different for CH 3 OH. Furthermore, the CH 3 OH temperatures do 
not show the same wide range of temperatures as the H 2 CO, even 


though these molecules have been shown to co-exist on large 
scales towards IRS7B (Lindberg & Jprgensen 2012). However, 
CH 3 OH observed in low-density outflow regions has been found 
to be extremely sub-thermally excited (Bachiller et al. 1995), and 
the CH 3 OH rotational temperature is not very sensitive to local 
increases in kinetic temperature (Bachiller et al. 1998). As seen 
in the CH 3 OH rotational diagram of the IRS7B survey (Fig. 4), 
this molecule shows a large scatter in the rotational diagram 
due to non-LTE effects. This shows that it is difficult to inter¬ 
pret CH 3 OH rotational diagrams containing as few as three data 
points. We therefore choose not to discuss the CH 3 OH rotational 
temperatures any further. 

The large spread in the H 2 CO temperatures might have been 
caused by optical depth effects, in particular for the H 2 CO 3 o 3 ^ 
2 o 2 line. We investigated this using the H^^CO 3 12 ^ 2ii tran¬ 
sition, which is detected towards seven of the studied sources. 
We assumed the ISM value for the ratio (69 + 6 ; Wilson 
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Fig. 7 . CH 3 OH rotational diagrams of the sources in the CrA survey where at least two CH 3 OH lines were detected and a positive rotational 
temperature could be fitted. The A lines are filled in the top half and the E lines are filled in the bottom half. Red square data points are 1.3 mm 
lines and blue circle data points are 0.9 mm lines, all observed with APEX. (A colour version of this plot is available in the online journal.) 


1999), an ortho-to-para ratio of 1.6, and excitation temperatures 
derived from the respective H 2 CO rotational diagrams. For a 
majority of the sources, we find that this line ratio corresponds 
to optical depths r of the H 2 CO Sqs ^ 2 o 2 line between 0.3 and 
0.8. For IRS2 and the SMM 2 on-source component the num¬ 
bers are higher ( 1.1 and 2 . 2 ), but it should be noted that these 
numbers are uncertain since the H^^CO 3 i 2 ^ 2 ii line is at a 
low S/N-level in these sources. The SMM 2 on-source value 
may also have errors originating in the separation between the 
on-source and outflow components. We conclude that the esti¬ 
mated H 2 CO optical depths are so low that the detected spread 
in temperatures must be real. 


The C-C 3 H 2 rotational temperatures are found to be constant 
at around 10-15 K for all sources. Non-LTE analysis using 
RADEX (van der Tak et al. 2007) shows that the observed line 
strengths are consistent with LTE and low optical depths, and 
we therefore expect that the measured rotational temperatures 
reflect the kinetic temperature of the C-C 3 H 2 gas. The H 2 CO 
temperatures are found to be strongly elevated compared to the 
C-C 3 H 2 temperatures (r(H 2 CO) ~ 30-50 K) in the sources close 
to R CrA (see Eig. 9). This is further discussed in Sect. 5. 


4.3. Molecular abundances in IRS7B and other protostars 
4.3.1. Rotational diagrams 

The only complex molecules detected towards IRS7B are 
CH 3 OH, CH 3 CCH, CH 3 CN, and CH 3 CHO, the latter two only 
tentatively detected. In typical hot corino sources, complex or¬ 
ganic molecules are present at comparably high abundance lev¬ 
els. To investigate whether the relatively low number of com¬ 
plex organic molecules detected towards IRS7B is significant, 
we compared relative abundances (and upper limits for non¬ 
detections) of organic molecules in IRS7B with the typical hot 
corino source IRAS 16293-2422, the WCCC source L1527, the 
quiescent source Bl-b, and the LI 157 outflow. To account for 
the uncertainty and difference in the sizes of the emission re¬ 
gions of the different sources, we normalised the abundances 
with that of CH 3 OH. We estimated the relative abundances by 
calculating the ratio of column densities as calculated from the 
rotational diagrams (see Table 3), thus assuming LTE and opti¬ 
cally thin lines. The large uncertainties on some of the column 
density measurements will, however, introduce large uncertain¬ 
ties in these ratios. 

To estimate upper limits on abundances of molecules not de¬ 
tected in IRS7B, we used Splatalogue and the line survey of 
IRAS 16293-2422 (Caux et al. 2011) to identify the strongest 
lines in the covered spectral bands. Assuming a temperature 
of 30 K and a line width of 2 km s“\ we calculated the abun¬ 
dance that would generate a line that would just reach the 3cr 
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Fig. 8 . C-C 3 H 2 rotational diagrams of the sources in the CrA survey where at least two C-C 3 H 2 lines were detected. The ortho lines are filled in 
the left half, the para lines are filled in the right half, and the lines that are blends of ortho transitions and para transitions are completely filled. All 
data points are 1.3 mm lines observed with APEX. (A colour version of this plot is available in the online journal.) 


detection limit in the APEX or ASTE observations of the band 
covering that certain line. In some ambiguous cases, several 
lines with similar line characteristics were tested, and the one 
generating the lowest upper limit was used. The lines used 
were for CH 3 OCHO the 227.560-227.564 GHz octuplet and for 
HCOOH the 220.038 GHz line. Eor CH 3 OCH 3 we found that 
the 225.599 GHz quadruplet produced the best upper limit in the 
APEX survey, but the 358.447-358.457 GHz octuplet covered 
by the ASTE survey (Watanabe et al. 2012) provided an upper 
limit three times lower, so this line was used instead. It is difficult 
to calculate upper limits of the abundances of complex organic 
molecules relative to CH 3 OH in the CrA source survey because 
the CH 3 OH excitation is poorly understood in these sources (see 
Sect. 4.2). 

The resulting abundance ratios or upper limits of a selec¬ 
tion of molecules measured towards IRS7B, the hot corino 
IRAS 16293-2422, the WCCC source L1527, the quiescent pro¬ 
tostar Bl-b, and the LI 157 outflow are shown in Table 6 and in 
Eig. 10. 

We note that several of the organic species are significantly 
less abundant towards IRS7B compared to IRAS 16293-2422, 
including CH 3 CN, CH 3 OCH 3 , CH 3 OCHO, and HCOOH, which 
makes the chemical parameters of IRS7B more similar to the 
LI 157 molecular outflow and the quiescent source Bl-b, al¬ 
though the CH 3 OH gas is much warmer towards IRS7B (28 + 1 K 
towards IRS7B, 10 ± 5 K towards Bl-b, and 12 ± 2 K to¬ 
wards the LI 157 outflow). Simpler organic species like H 2 CO, 


H 2 CCO, and HNCO have similar abundances towards IRS7B 
and IRAS 16293-2422, whereas unsaturated hydrocarbons and 
cyanides (C-C 3 H 2 , CN, HC 3 N, and C 2 H) are over-abundant to¬ 
wards IRS7B compared to IRAS 16293-2422. The unsaturated 
hydrocarbons C-C 3 H 2 and C 2 H are, however, under-abundant to¬ 
wards IRS7B compared to the WCCC source L1527, and C 4 H 
is non-detected towards IRS7B (Sakai et al. 2009a). This shows 
that IRS7B is not a WCCC source. The given upper limits, as 
well as the values calculated for H 2 CCO, HNCO, CH 3 CN, and 
CH 3 CHO, were calculated under the assumption that the tem¬ 
perature of these species is 30 K. A lower temperature would 
increase the ratios for these species. 

Interferometric observations of CH 3 OH show that this emis¬ 
sion is present on large scales around IRS7B, in particular in 
two extended ridges north and south of IRS7B (Lindberg & J 0 r- 
gensen 2012), which are partially covered by the APEX and 
ASTE beams. These beam-averaged abundance ratios therefore 
probably reflect the composition of that large-scale molecular 
gas heated by R CrA (Lindberg & Jprgensen 2012), rather than 
the compact inner envelope of IRS7B. 

4.3.2. Radiative transfer modelling 

We used the radiative transfer code RATRAN (Hogerheijde 
& van der Tak 2000) to model the line strengths of some 
of the molecular species detected towards IRS7B. The one¬ 
dimensional envelope and temperature model made by Lindberg 
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Table 4 . Rotational diagram parameters for H 2 CO in the survey of CrA sourees. Sources where the fit could not be performed are excluded. The 
H 2 densities n and column densities N were measured by the method described in Appendix F, which can only be made for the sources where both 
the 3 o 3 ^ 2 o 2 and the 5 o 5 ^ 4 o 4 lines have been detected. To facilitate comparison between IRS7B and the other sources, only the lines covered 
by the source survey are included, which explains why IRS7B shows different values here and in Table 3. 


Source name 

Tj-3^ 

[K] 


b 

[K] 

7’av‘= 

[K] 


n 

[cm“^] 


[cm“^] 

IRS7B 

48.8 + 

6.2 

43.8+ 2.7 

45.3 + 

3.4 

inx 

10^ 

1.0 X 

Haas4 

< 24.3 








IRS2 

18.9 + 

1.3 


18.9 + 

1.3 




IRS5A 

29.8 + 

2.4 


29.8 + 

2.4 

4.0 X 

10^ 

3.7 X 

IRS5N 

28.4 + 

2.1 


28.4 + 

2.1 

3.9 X 

10^ 

3.9 X 10^3 

IRSl 

29.5 + 

2.4 


29.5 + 

2.4 

4.6 X 

10^ 

2.8 X 1013 

IRS7A 



43.6+ 3.1 

43.6 + 

3.1 




CrA-24 source 

40.3 + 

4.3 

28.0+ 3.1 

33.2 + 

2.7 

1.3 X 


2.6 X 1013 

CrA-24 outflow 

68.6 + 

12.9 

103.7 + 68.1 

74.2 + 

34.7 

e 


e 

SMM 2 source 

40.4 + 

5.3 


40.4 + 

5.3 

6.3 X 

10^ 

6.6 X 1013 

SMM 2 outflow 

52.6 + 

7.3 


52.6 + 

7.3 

<2 X 

10^^ 

>6 X 1013^ 

CX042 

< 87.5 








CrA-44 

< 21.2 








CrA-33 

< 54.5 








VV CrA 

32.7 + 

3.7 


32.7 + 

3.7 





Notes. Rotational temperatures of the transitions with /u = 3. Rotational temperatures of the transitions with /u = 5. Weighted average 
of the two previous, if they are both measured, otherwise the same as the only one measured. The total H 2 CO column density assuming an 
ortho-to-para ratio of 1.6. For these sources, the column density and H 2 density could not be constrained from the 3 o 3 ^ 2 o 2 / 5 o 5 ^ 4 o 4 ratio. 

Table 5 . Rotational diagram parameters for CH 3 OH and C-C 3 H 2 in the survey of CrA sources. Sources where neither of the fits could be performed 
are excluded. To facilitate comparison between IRS7B and the other sources, only the lines covered by the source survey are included, which 
explains why IRS7B shows different values here and in Table 3. 


Source name 

T 

[K] 

CH 3 OH 

[cm- 2 ] 


T 

[K] 

C-C 3 H 2 

ATb 

[cm-3] 

IRS7B 

13.5+ 0.6 

(4.5 + 1.0) X 

1014 

16.5 + 2.6 

( 8.6 + 3.2) X 1013 

Haas4 




10.6+1.7 

(1.1 +0.5) X 1013 

IRS2 

< 16.2 

> 1.2 X 

1013 



IRS5A 

18.1+ 2.6 

(4.5 +2.1) X 

1013 

9.9+ 0.8 

(3.7 ± 0.9) X 1013 

IRS5N 

17.7+ 2.2 

(5.3 + 2.2) X 

1013 

10.9 + 0.9 

(3.3 + 0.8) X 1013 

IRSl 

18.9+ 2.4 

(6.7 + 2.6) X 

1013 

8.9 + 1.0 

(1.3 + 0.5) X 1013 

IRS7A 

12.2+ 0.9 

(6.4 + 3.0) X 

1014 



CrA-24 




13.0+1.7 

(2.1 ±0.7) X 1013 

CrA-24 source 






CrA-24 outflow 

46.3 + 20.2 

(4.4 + 2.7) X 

1013 



SMM 2 




14.7 + 2.9 

(1.2 ± 0.5) X 1013 

SMM 2 source 

30.6+11.4 

(1.4+ 1.0) X 

1013 



SMM 2 outflow 

19.9+ 2.3 

(1.6 + 0.6) X 

1014 



CrA-44 




11.4 + 0.9 

(6.0+ 1.4) X 1013 

CrA-33 

< 14.4 

> 2.2 X 

1013 

< 10.7 

>5.8 xlO" 

VV CrA 




12.5 + 1.8 

(7.8±2.9)x lOii 


Notes. The CH3OH column density assuming an A/E ratio of 1. The C-C3H2 column density assuming an ortho-to-para ratio of 3. 


& J 0 rgensen (2012) and refined by Lindberg et al. (2014a) was 
employed. This is an envelope with a power-law density profile 
(p ~ heated both by a central object and from the outside 
by irradiation corresponding to that of the nearby Herbig Be star 
R CrA. We used the RATRAN models to calculate abundances 
for 13 important species for which most or all lines are expected 
to be optically thin (see Jprgensen et al. 2004). The model abun¬ 


dances were kept constant as a function of radius for all species 
except CN, since the modelled species are not expected to be 
heavily affected by hot corino emission. On the other hand, we 
let the abundance of the photo-destruction product CN take the 
shape of a step function, allowing for a higher CN abundance 
in the outer envelope, where the gas is more affected by photo¬ 
destruction from UV radiation originating in R CrA. 
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Table 6. Molecular abundances relative to CH 3 OH (NiXyNiCR^OH)) towards IRS7B and other protostellar sources as established from rotational 
diagrams. See also Fig. 10. 


Molecule 

IRS7B 

116293^ 

L1527'’ 

Bl-b core" 

LI 157 outflow‘d 

CH 3 OH 

= 1 

= 1 

= 1 

= 1 

= 1 

H 2 CO 

0.62 

0.20 

0.67 



H 2 CCO 

0.047 

0.041 




HNCO 

0.032 

0.039 

0.052 



CH 3 CHO 

0.10 

0.17 


0.012 

0.006 

HCOOH 

<0.03 

0.21 




CH 3 CN 

0.006 

0.033 




CH 3 OCH 3 

<0.03 

0.80 


< 0.008 


CH 3 OCHO 

<0.16 

1.3 

<0.08 

0.023 

0.018 

C-C 3 H 2 

0.053 

0.008 

0.21 



CH 3 CCH 

0.67 

0.87 

0.95 



CN 

1.1 

0.023 

> 0.76" 



HC 3 N 

^4 

0.006 

0.43 



C 2 H 

1.8 

0.057 

15.6 




Notes. IRAS 16293-2422 (van Dishoeck et al. 1995; Schoier et al. 2002; Cazaux et al. 2003). Sakai & Yamamoto (2013) and references 
therein. Observations were performed with several different single-dish telescopes with beam sizes between 10" and 30". Since the source is 
extended (Sakai et al. 2010), we used the beam-averaged column densities. Oberg et al. (2010). Bachiller & Perez Gutierrez (1997); Arce 
et al. (2008). Agiindez et al. (2008). 
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Fig. 10. Graphical representation of molecular abundances relative to the CH 3 OH abundance as given in Table 6 for IRS7B, the hot corino source 
IRAS 16293-2422, the WCCC source L1527, the quiescent protostar Bl-b, and the LI 157 outflow. See the notes of Table 6 for references. For 
a discussion on the classification of gas and grain formation of the organic species, see Charnley (1997, 2001). (A colour version of this plot is 
available in the online journal.) 


In addition to the APEX line data of IRS7B reported in 
this work, we used C-C 3 H 2 , C^^O, sO, CN, HC^^O+, and 
DCN line data from the ASTE survey of IRS7B (Watanabe et al. 
2012 ). 

Molecular data files from the LAMDA database (Schoier 
et al. 2005) were employed when available. In addition, the same 
molecular datafiles as were used by Jprgensen et al. (2004) were 
used for CN and HC 3 N, since data for hyperfine lines of CN 


and high-/ transitions of HC 3 N are not available in the LAMDA 
database. Eor isotopologue species where no LAMDA data files 
are available, the collisional data of the main isotopologue were 
used, but the line parameters of the isotopologue were acquired 
from the CDMS database (Muller et al. 2001). This was done 
for *3CS, c^^S, DCN, and DNC. 

The line widths and LSR velocities were estimated from 
Gaussian fits to the spectral lines. The abundance was left as 
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Fig. 9. Rotational temperatures for a) H 2 CO, b) CH 3 OH, and 
c) C-C 3 H 2 as functions of distance to R CrA. Upper limits are shown 
in grey. In the H 2 CO plot, CrA-24 and SMM 2 have two data points 
each, one for the on-source component (blue) and one for the outflow 
component (red), and the same applies to the CH 3 OH plot when the 
rotational temperature could be measured. The CH 3 OH rotational tem¬ 
perature most likely does not reflect a physical temperature due to non- 
LTE excitation in combination with few data points. (A colour version 
of this plot is available in the online journal.) 


a free parameter, and the best value was found through a least-;^^ 
fit. The error used for the calculations is the combination of 
the rms error and a calibration uncertainty of 10 % for both the 
APEX and the ASTE observations. 

Table 7 lists the resulting modelled abundances together with 
data for other low-mass embedded protostellar sources in the lit¬ 
erature. The abundance of the optically thick species CO, CS, 
and HCO^ were calculated from optically thin isotopologues us¬ 
ing isotope ratios from the literature: C, O, and isotope ratios 
were taken from Wilson (1999) and isotope ratios from Chin 
etal. (1996). 

A contour plot of the reduced-;^^ values as function of the 
CN abundances is shown in Eig. 11. The fit clearly is relatively 
insensitive to the inner abundance (Xinner ^ 10 “^), but fixes the 
outer abundance to Xouter ~ few x 10“^. 



Fig. 11. Reduced-;^^^ fit to the CN abundance. The solid lines show the 
Icr, 2cr, and 3 cr confidence levels. The dashed lines show the reduced- 
values in steps of 5. The border between inner and outer abundance 
lies where the temperature falls below 30 K, which occurs at 370 AU. 
The blue dotted line shows the values for equal inner and outer abun¬ 
dances. (A colour version of this plot is available in the online Journal.) 

While the IRS7B abundances of many of the more common 
molecules (CO, HCO^, SO, and CS) are comparable to those in 
other sources, we find that the CN-bearing species (CN, DCN, 
DNC, and HC 3 N) and C-C 3 H 2 are strongly enhanced in IRS7B 
compared to the other protostellar envelopes. This, in combina¬ 
tion with the low upper limits of complex organic molecules dis¬ 
cussed above, demonstrates that the chemical properties of the 
gas around IRS7B strongly differs from those of most other pro¬ 
tostars studied before. Previous studies have shown that small 
hydrocarbon molecules such as C-C 3 H 2 and C 2 H are enhanced 
in PDR regions. In high-UV flux regions this is a result of gas- 
phase chemistry (Cuadrado et al. 2015), but in regions with more 
moderate UV flux levels it is a consequence of photodestruction 
of PAHs (Montillaud et al. 2013). For further references and 
discussion, see also Guzman et al. (2015). 

4.4. Isotopic fractionation in CrA 

The [C^^O]/[C^^O] abundance ratio in IRS7B is found to be 5.6, 
considerably higher than the expected value in the local ISM, 
3.6 (Wilson 1999). It should, however, be noted that the C^^O 
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Table 7 . Results of RATRAN models of fractional abundances of certain molecules in R CrA IRS7B and similar modelling results from the 
literature for other protostellar sources. The reduced-;^^^ values for IRS7B assume a 10% calibration error and can only be computed for species 
with at least two spectral lines. 


Molecule 

IRS7B 

Abundance 

xid 

116293^ 

Abundance 

Class OO 
Abundance 

Class lO 
Abundance 

Prestellar^ 

Abundance 

CO" 

1.4 X 10-4 


3.3 X 

10-5 

2.1 X 10-5 

1.1 X 10-4 

1.4 X 10-5 

c'^o 

5.7 X 10-^ 

6.9 

1.6 X 

10-^ 




C'^0 

3.2 X 10-’ 


6.2 X 

10-^ 




HCO+" 

3.8 X 10-® 


1.4 X 

10-1 

2.2 X 10-1 

1.1 X 10-* 

8.0 X 10-10 

n 

00 

0 

+ 

6.8 X 10-12 

3.6 

6.4 X 

10-12 




so 

1.9 X lO-ii 

33.6 

4.4 X 

10-1 

2.0 X 10-1 

2.3 X 10-1 

1.5 X 10-1 

cs^ 

6.7 X 10-11 


3.0 X 

10-1 

1.5 X 10-1 

4.3 X 10-1 

1.8 X 10-1 

‘^cs 

8.8 X 10-" 







C34S 

2.1 X 10-10 

16.2 

1.2 X 

10-10 




c^^s 

6.8 X 10-11 







CN outer^ 
CN inner^ 

4.5 X lO-o 
< 1 X 10-1 

8.5 

8.0 X 

10-11 

6.9 X 10-10 

3.0 X 10-1 

4.8 X 10-1 

DCN 

7.2 X 10-11 

14.1 

1.3 X 

10-11 

1.4 X 10-11 


7.7 X 10-12 

DNC 

6.8 X 10-11 


4.2 X 

10-12 




HNCO 

1.2 X 10-10 

0.002 

1.3 X 

10-10 




HC 3 N 

3.3 X 10-1 

11.2 

1.5 X 

10-10 

3.5 X 10-10 

3.1 X 10-1 

8.2 X 10-10 

C-C3H2 

3.0 X 10-10 

12.9 

1.6 X 

10-11 





Notes. IRAS 16293-2422, Schoier et al. (2002). Average of sources in Jprgensen et al. (2004). The abundances for these optically thick 
species have been calculated from the optically thin isotopologues assuming local ISM isotope ratios (Wilson 1999; Chin et al. 1996). The CN 
fit was performed with a step in the abundance profile at 370 AU (where T ^ 30 K). 


fit includes both the 3-2 and 2-1 lines, whereas the fit 
only includes the 2-1 line. The uncertainty of the large-scale 
CO distribution and the differing beam sizes could contribute 
to the unusual isotope ratio, and if the ratio is calculated only 
with the 2-1 lines, a [C^^O]/[C^^O] ratio of 4.0 is reached. The 
[^"^CS]/[C^^S] ratio 2.4 is slightly lower than the local ISM value 
2.9, and the [^'^CS]/[^^CS] ratio 3.1 is lower than the local ISM 
value 6.3 + 1.0. 

We also investigated the D/H ratio of H 2 CO towards IRS7B. 
As mentioned in Sect. 4.1, the excitation conditions of the H 2 CO 
isotopologues is somewhat uncertain. Since the D 2 CO rotational 
diagram fit is much better than the H^^CO fit, we applied the 
D 2 CO rotational temperature to the H^^CO excitation diagram 
to acquire a better estimate on the H^^CO column density. As¬ 
suming the local ISM value for the ratio (69 + 6 , Wilson 

1999), we estimated the [D 2 CO]/[H 2 CO] ratio towards IRS7B 
to 0.018 + 0.003. The calculated ratio is within the errors of the 
value of Watanabe et al. (2012), who calculated the ratio for tem¬ 
peratures fixed at 15 K, 20 K, and 25 K. Their [D 2 CO]/[H 2 CO] 
ratio for 25 K is 0.021 + 0.003. They also observed HDCO, 
and calculated the [HDCO]/[H 2 CO] ratio to 0.052 + 0.008 as¬ 
suming the same temperature. We compared these values with 
[D 2 C 0 ]/[H 2 C 0 ] and [HDCO]/[H 2 CO] ratios in other sources. 
Although the [D 2 CO]/[H 2 CO] value is elevated in comparison 
with what was found in Orion KL ([D 2 CO]/[H 2 CO] = 0.003; 
Turner 1990), the level of deuteration is lower than what is 
found in typical hot corino sources and other Class 0 proto¬ 
stars (Parise et al. 2006; e.g. [D 2 CO]/[H 2 CO] = 0.05 + 0.008 
in IRAS 16293-2422). This could suggest that grain-surface 
chemistry is less important for the formation of H 2 CO in IRS7B 
than in other low-mass protostellar sources, but the lower level 
of deuteration could also be an effect of the higher large-scale 
temperature due to the external irradiation. However, while the 


[D 2 C 0 ]/[H 2 C 0 ] and [HDCO]/[H 2 CO] ratios are biased by the 
D/H ratio in the precursor molecules, the incremental fraction¬ 
ation ratio [HDC 0 ]/[H 2 C 0 ]:[D 2 C 0 ]/[HDC 0 ] only depends on 
the formation mechanisms for H 2 CO (Turner 1990; Rodgers & 
Chamley 2002). The incremental fractionation ratio of H 2 CO in 
IRS7B is 0.15 + 0.04, assuming the [D 2 CO]/[H 2 CO] from this 
work and the [HDCO]/[H 2 CO] ratio of Watanabe et al. (2012) 
at 25 K. This is a low value, but still within the errors of the hot 
corino values (Parise et al. 2006), so we cannot draw any clear 
conclusions as to whether there is a difference in the formation 
path of H 2 CO in IRS7B versus hot corinos or not. 

Watanabe et al. (2012) also measured the DCN/HCN ratio, 
which was found to be 0.010+0.001, similar to the IRAS 16293- 
2422 value (van Dishoeck et al. 1995). 

We detect and towards several of the sources in 
CrA. In most sources, the ^^CN lines are only observed in the 
3^2 rotational transitions, while the isotopologue species are 
observed in the 2^1 rotational transitions. The uncertainty in 
the excitation conditions makes the estimation of isotope ratios 
difficult. For IRS7B, we have 2^1 observations also for the 
^^CN lines, and (excluding the strongest ^^CN lines, which are 
expected to be optically thick) we find a CN/C^^N value of 229 ± 
50, which is within the errors of the ratios calculated 

from IRS7B observations using H^^CN/HC^^N ratios (287 ± 36) 
and HN^^C/H^^NC ratios (259 ± 34) by Wampfler et al. (2014). 

Faint lines are detected towards IRS7A (but not 

IRS7B). The resulting [SO]/[34 sq] to be 11 ± 10, 

which is almost within the errors of the ISM value ~ 22 (Wilson 
1999). 
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Fig. 12. 1-D radiative transfer model of heating from R CrA (dashed 

line) plotted together with the H 2 CO rotational temperatures measured 
in the CrA sources as a function of distance to R CrA. For CrA-24 and 
SMM 2, only the on-source components are shown. Upper limits are 
shown in grey. See also Fig. 9a. 



5. Discussion 

In Sect. 4.2, the H 2 CO rotational temperatures in the survey of 
CrA sources were found to be highest close to the Herbig Be 
star R CrA (Fig. 9a). This suggests that the sphere of thermal 
influence of R CrA (see Adams & Shu 1985) stretches out to 
radii greater than 10000 AU. To study if this is plausible, ra¬ 
diative transfer modelling was employed. Using the dust radia¬ 
tive transfer code Transphere (Dullemond et al. 2002), we con¬ 
structed a one-dimensional model of the molecular cloud around 
R CrA. The density was assumed to be constant at ^ = 10^ cm“^ 
throughout the cloud, starting at a distance of 200 AU from 
R CrA. The luminosity of R CrA was set to 144 L© (see Lind- 
berg et al. 2014b for a thorough discussion of the luminosity of 
R CrA). For this configuration, Transphere calculates the tem¬ 
perature in the region self-consistently as a function of radius. 
The model is shown compared to the data in Fig. 12, and we 
conclude that it provides a good fit to the H 2 CO rotational tem¬ 
perature data points within a radius of ~ 30 000 AU from R CrA. 
We attribute the deviations outside the error bars to irregularities 
in the distribution of gas in the molecular cloud, with the excep¬ 
tion for the source VV CrA. 

The protostar VV CrA has a projected distance of more than 
0.5 pc (100000 AU) to R CrA, but still has a comparably high 
H 2 CO rotational temperature (33 + 4 K). This cannot be ex¬ 
plained by the irradiation from R CrA. However, VV CrA is 
a fairly luminous protostar, with a higher flux at 5.8 [im than 
R CrA (Peterson et al. 2011). Using 2MASS, Spitzer, IRAS, 
PACS, and SCUBA data of VV CrA (Sicilia-Aguilar et al. 2013, 
and references therein), we estimate the bolometric luminosity 
of the source to 13L0, which is not sufficient to explain the mea¬ 
sured H 2 CO temperature. However, by the use of the Robitaille 
database of SED models for YSOs (Robitaille et al. 2006) and 
their online^ fitting tool (Robitaille et al. 2007), we find that the 
best matches to the SED are protostellar models with luminosi¬ 
ties > 4 OL 0 . Such a luminosity can heat the gas to 30 K on scales 
of several 1000 AU, and the elevated temperature in the envelope 
of VV CrA could thus be attributed to internal irradiation from 


^ http://caravan.astro.wise.edu/protostars/ 


VV CrA itself. It should, however, be noted that the SED mod¬ 
els of Robitaille et al. (2006) work best with protostars at later 
stages of evolution. 

The C-C 3 H 2 rotational temperatures measured in the source 
survey are much lower than the H 2 CO temperatures and inde¬ 
pendent of the distance to R CrA. This can be explained by the 
C-C 3 H 2 gas being more compact than the H 2 CO gas, and thus 
more shielded against external irradiation. SMA interferometry 
of the region around IRS7B shows C-C 3 H 2 emission that is more 
concentrated towards the low-mass YSOs than the H 2 CO emis¬ 
sion is, which supports this hypothesis (Lindberg & Jprgensen 
2012 ), but H 2 CO and C-C 3 H 2 mapping of additional sources is 
required to confirm this explanation. 

We also investigated the chemical composition of the molec¬ 
ular gas detected towards IRS7B. A simple comparison of spec¬ 
tra from three different star-forming regions - a hot core, a hot 
corino, and IRS7B (Eig. 3) - suggests that the chemistry of the 
low-mass hot corino source IRAS 16293-2422 has more in com¬ 
mon with the high-mass hot core G327.3-0.6 than with the low- 
mass protostar IRS7B. Towards IRS7B, the PDR-tracing CN 
lines dominate the spectrum completely, whereas the two other 
sources have spectra dominated by various complex organic 
molecules. The only complex molecules detected in IRS7B are 
CH 3 OH, CH 3 CCH, CH 3 CN, and CH 3 CHO (the latter two only 
tentatively detected). CH 3 OCH 3 and CH 3 CN are non-detected 
also in high-resolution ALMA observations of the inner enve¬ 
lope of IRS7B (Lindberg et al. 2014a). We also find that the rela¬ 
tive abundances of several complex organic molecules in IRS7B 
are significantly lower than what is found in typical hot corino 
sources. However, simple cyanides (CN and HCN) and the un¬ 
saturated hydrocarbons C-C 3 H 2 and HC 3 N are over-abundant in 
IRS7B. This makes IRS7B similar to the warm carbon-chain 
chemistry sources, but considering that the C 2 H abundance is 
an order of magnitude higher in LI527 than in IRS7B and that a 
very low upper limit of the C 4 H abundance in IRS7B was mea¬ 
sured by Sakai et al. (2009a), IRS7B is not likely to be a typical 
WCCC source. Rather, the detected molecular abundances re¬ 
flect a PDR-like chemistry resulting from the strong irradiation 
onto the molecular envelope from R CrA. 

The CH 3 OH abundance observed towards IRS7B is only 
around 3-6 x 10“^ (we calculated this indirectly using the ro¬ 
tational diagram column density of CH 3 OH and the rotational 
diagram column densities and RATRAN abundances of HNCO 
and C-C 3 H 2 ). This is very low when compared to hot corino 
abundances (3 x 10“^ in IRAS 16293-2422; Cazaux et al. 2003) 
and could be explained by the elevated temperature in the region, 
which inhibits grain formation of CH 3 OH due to evaporation 
of both CO and H atoms. The low abundances of CH 3 OCHO 
and CH 3 OCH 3 relative to CH 3 OH can be explained if we as¬ 
sume that these molecules are generally formed through gas- 
phase reactions (see, e.g., Charnley 1997; Peeters et al. 2006). 
In this scenario, a CH 3 OH abundance < 10“^ is too low for an 
efficient formation of these molecules (Taquet et al., in prep.). 
H 2 CCO and CH 3 CHO are, on the other hand, more likely to 
form on grains (Charnley 2001) and should form at a rate sim¬ 
ilar to CH 3 OH, which could explain why the H 2 CCO/CH 3 OH 
and CH 3 CHO/CH 3 OH ratios are similar to what is found in 
IRAS 16293-2422. Another possibility is that the tentatively 
detected complex molecules could exist in a more shielded in¬ 
ner region, at temperatures at the same levels as those of the 
hydrocarbon molecules (10-15 K), similar to the properties of 
the complex organic molecules detected in the prestellar core 
L1689B (Bacmann et al. 2012). If, on the other hand, we assume 
a grain-formation scenario for all complex organic molecules. 
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Fig. 13. As a proxy for the CH 3 OH/H 2 CO abundance ratio, we here 
show the line ratio between the CH 3 OH 42 ^ 3i, E line at 218.440 GHz 
and the H 2 CO 3 o 3 ^ 2 o 2 line at 218.222 GHz in the CrA sources as 
a function of the H 2 CO rotational temperature. CrA-24 and SMM 2 
are shown with two data points each, one for the on-source component 
(blue) and one for the outflow component (red). (A colour version of 
this plot is available in the online journal.) 


the elevated temperatures should enhance the formation of com¬ 
plex organic molecules as long as precursor molecules such as 
H 2 CO exist on the icy grains (Garrod & Herbst 2006). The 
CH 3 CN found towards protostellar envelopes is believed to 
be formed from gas-phase reactions between HCN and CH 3 
(Rodgers & Chamley 2001; Wang et al. 2010). The low CH 3 CN 
abundance and high CN and DCN abundances observed towards 
IRS7B are thus puzzling, but might be due to a relative lack of 
CH 4 , the precursor species of CHJ. The exact reason for a low 
CH 3 CN abundance needs to be further investigated. 

Laboratory experiments (Watanabe et al. 2004; Fuchs et al. 
2009) and numerical simulations (Cuppen et al. 2009) have sug¬ 
gested that the hydrogenation reactions of CO forming H 2 CO 
and CH 3 OH on CO-H 2 O ices are strongly temperature depen¬ 
dent, with the CH 3 OH/H 2 CO ratio increasing significantly at 
relatively small increases in temperature. Since the CH 3 OH ro¬ 
tational diagrams of the source survey have too few data points 
to be reliable (see Sect. 4.2), we did not use the CH 3 OH/H 2 CO 
abundance ratio, but instead used the ratio between the CH 3 OH 
42 ^ 3i, E line at 218.440 GHz and the H 2 CO 3 o 3 ^ 2 o 2 line 
at 218.222 GHz as a more useful proxy for the CH 3 OH/H 2 CO 
abundance ratio. We compare this ratio with the H 2 CO rota¬ 
tional temperature in Fig. 13 and see a positive trend between 
the CH 3 OH/H 2 CO ratio and the H 2 CO temperature. We sug¬ 
gest that the high level of irradiation has caused most of the CO 
to evaporate from the ices in the whole region around R CrA, 
which has caused the relatively low CH 3 OH abundances (see 
above). However, the low amounts of CO that are still left on the 
grain surfaces will form CH 3 OH more efficiently in the warmest 
protostellar envelopes, causing the trend in Fig. 13. This effect 
should be strongest at T < 20 K (Watanabe et al. 2003), but the 
trend observed above 20 K could be explained if the CH 3 OH 
now observed in the gas phase was formed on the icy grains at a 
time when R CrA was less luminous than today. 

Lindberg & Jprgensen (2012) found that the embedded pro¬ 
tostellar sources in the R CrA cloud are encompassed by dense 
molecular gas traced by ridges of H 2 CO on scales of at least 
5 000 AU. A large part of the emission observed in the APEX 


beam centred at IRS7B is expected to originate in these ridges. 
Towards CrA-24 and SMM 2, we find signs of large-scale H 2 CO 
and CH 3 OH gas at velocities somewhat lower than the typi¬ 
cal vlsr of the sources in the CrA region. This is consistent 
with CH 3 OH emission extending towards the northeast from 
the R CrA cloud observed in SEST CH 3 OH maps (Miettinen 
et al., in prep.). A schematic drawing of the distribution of dense 
molecular gas, dust emission, and YSOs in the vicinity of R CrA 
is shown in Fig. 14. A direct comparison with the typical hot 
corino sources might therefore be argued against, because in 
such sources the emission from complex organic molecules orig¬ 
inates in relatively compact regions close to the central source 
(e.g. Kuan et al. 2004; Bottinelli et al. 2004b, 2008; Bisschop 
et al. 2008). However, Lindberg et al. (2014a) separated the faint 
CH 3 OH emission originating in the IRS7B central source from 
the surrounding emission using high-resolution ALMA interfer¬ 
ometry data. The faint on-source CH 3 OH lines can either be ex¬ 
plained by a very low CH 3 OH abundance in the inner envelope 
(X(CH 30 H) ~ 10 “^^; an order of magnitude lower than the al¬ 
ready low large-scale CH 3 OH abundance), or by a more typical 
CH 3 OH abundance in the inner envelope (X(CH 30 H) ~ 10“^) in 
combination with a flattened density profile in the inner 100 AU 
caused by the presence of a Keplerian disc. This latter case 
would still allow for a higher CH 3 OH abundance and the pres¬ 
ence of more complex organic molecules in the inner envelope of 
IRS7B at similar abundances as found in hot corinos. Chemical 
modelling of the irradiated envelope and ALMA observations 
with higher sensitivity and resolution will help clarify the causes 
of the chemistry observed towards IRS7B and the inner structure 
of the source. 


6. Conclusions 

We have performed an unbiased line survey of the Class O/I 
source R CrA IRS7B to investigate its chemical and physical 
properties. We also presented a survey of the H 2 CO, CH 3 OH, 
and C-C 3 H 2 emission for the full sample of embedded YSOs 
in the Corona Australis region. The observations were carried 
out with the (sub)millimetre single-dish APEX telescope, but we 
also used supplementary data from the ASTE telescope. These 
are our main conclusions: 

1. In comparison with other well-studied Class O/I sources, the 
deeply embedded YSO IRS7B shows a different chemistry. 
In particular, very strong CN emission is seen, which in¬ 
dicates a PDR-like chemistry. The only signs of complex 
organic molecules or large carbon-chain species detected 
towards IRS7B are CH 3 OH, CH 3 CCH, and faint CH 3 CN 
and CH 3 CHO emission. The CH 3 OH abundance is also 
very low (3-6 x 10“^) in comparison with to other embed¬ 
ded protostars. The values or estimated upper limits of the 
abundances of several complex organic species relative to 
CH 3 OH are much lower than the relative abundances de¬ 
tected towards typical hot corino sources, while species such 
as CN, HC 3 N, and C-C 3 H 2 show strongly elevated relative 
abundances. High-resolution ALMA observations have also 
failed to detect complex organics on smaller scales (Lindberg 
et al. 2014a). 

2. The molecules detected towards R CrA IRS7B can be 
grouped into two categories by their rotational temperatures: 
the warm carbon-chain molecules and nitrogen-bearing 
species show temperatures around 10-17 K, whereas the 
organic molecules and sulphur-bearing species are warmer 
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Fig. 14. Schematic drawing of the R CrA cloud and its surroundings. Protostars are shown in yellow. Red structures indicate extended H 2 CO 
and CH 3 OH emission, and blue striped structures show extended emission from hydrocarbon and nitrogen-bearing species (Miettinen et al., in 
prep.; Lindberg & Jprgensen 2012). The brown dashed contour shows extended high-/ (/ ~ 20) CO emission, and the black dotted contours show 
extended 110 |im dust continuum emission not associated with the point sources (from Herschel PACS observations; Lindberg et al. 2014b). The 
solid black circle shows the APEX beam at 218 GHz centred at IRS7B, and the dashed black circle shows the APEX beam at 363 GHz. Large 
parts of the region have not been covered by all observation modes, and thus the map is probably incomplete. This map covers only a small part 
of the sources in the CrA source survey (cf. Eig. 1). (A colour version of this plot is available in the online journal.) 


at 19-40 K. Molecules from the two groups have pre¬ 
viously been found to have distinct spatial distributions 
in SMA/APEX interferometry observations of the region 
(Lindberg & Jprgensen 2012). 

3. The H 2 CO rotational temperatures of the sources in the CrA 
survey near the Herbig Be star R CrA are elevated to 30- 
50 K, and the H 2 CO temperature is also found to decrease 
with the distance to R CrA. This temperature gradient is 
consistent with a 1-D radiative transfer model with R CrA 
as the heating source. The H 2 CO temperatures of the out¬ 
flow components of the sources with two velocity compo¬ 
nents are higher than the corresponding on-source compo¬ 
nent temperatures. The C-C 3 H 2 rotational temperatures are 
on the other hand found to be fairly constant around 10-15 K 
for all sources where it could be measured. This could be a 
shielding effect as a result of the C-C 3 H 2 being present at 
much smaller spatial scales than the H 2 CO. 

4. The CH 3 OH/H 2 CO is found to increase with the H 2 CO ro¬ 
tational temperature. This can be explained by the higher 
temperature increasing the efficiency of hydrogenation reac¬ 
tions of CO, as suggested by laboratory studies. 

We have shown that the chemistry of protostellar envelopes 
can be strongly influenced by external irradiation. Follow-up 
studies with interferometric observations are needed to reveal 
the different origins of various organic molecules and unsatu¬ 
rated hydrocarbons in protostellar envelopes. To understand why 
the abundance levels of some but not all organic molecules are 
lowered, their exact formation paths need to be studied using 
astrochemical modelling, in particular assuming elevated tem¬ 
peratures on large scales. 
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Appendix A: Baseline subtraction 

As noted in Sect. 2, the spectra observed by APEX exhibit unstable quasi-sinusoidal baselines (see Vassilev et al. 2008). When 
investigating the spectra, we found that polynomial and/or sinusoidal baselines were not giving desirable results. To account for 
this, we developed a relatively advanced baseline-fitting algorithm, computing running-mean baselines, which is described here. 

In principle, the fitted baseline is a boxcar smoothing of the spectrum (where line channels have been removed) with weights 
increasing towards the centre of the box: For each channel in the spectrum, the algorithm calculates the weighted mean of the 
line-free channels in a box around the channel. The weights have a Gaussian distribution around the central channel of the box. 
The width of this Gaussian must be adjusted so that all irregularities in the baseline are removed but no real spectral features are 
removed. In addition to the quasi-sinusoidal baselines, strong atmospheric lines arising from the different elevations of the position¬ 
switching on- and off-positions made it necessary to introduce discontinuities in the fitted baseline. By this method, a piecewise 
smooth but in general non-analytical baseline can be fitted to the data (see examples in Fig. A.l). Some artefacts of the strongest of 
these atmospheric lines still remain after the baseline subtraction and appear as broad spectral line features in the resulting spectra. 
In this dataset, they are easy to distinguish from the source lines due to the large difference in line widths. 
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Frequency [GHz] Frequency [GHz] 


Fig. A.l. Examples of spectral data before they were treated by the running-mean baseline-fitting routine (blue) and the fitted baseline that was 
subsequently subtracted (red). The right spectrum has a discontinuity in the fitted baseline that is due to an atmospheric line. (A colour version of 
this plot is available in the online journal.) 


Appendix B: Observed spectral line parameters 

The integrated intensities and other parameters of the detected spectral lines are listed in the tables of this appendix. The values of 
vlsr, Av, and Tpeak are results of Gaussian fits. In the cases where no values are given, Gaussian fits could not be performed to that 
line because of irregular line shapes. 

Table B.l. Observed line strengths towards IRS7B, both for the 1.3 mm unbiased 
line survey and for the IRS7B part of the source survey, which also covered the 
0.9 mm and 0.8 mm windows. 


Frequency^ 

[MHz] 

Molecule 


Transition 


VlSR 
[km s"^] 

Av 

[km s"^] 

7’mb 

[K] 

j 7'mbrfv 
[Kkms-*] 

217238.538 

DCN 


7 = 3^2 


6.1 

1.8 

1.04 

2.01 + 0.07 

217467.150 


N 

= 2^ l,/ = 5/2^3/2,Fi = 

3^2,F = 4^3 

5.8 

0.9 

0.24 

0.24 + 0.05 

217467.150 


N 

= 2^ l,/ = 5/2^3/2,Fi = 

3 ^2,F = 3 ^2 





217822.148 

C-C 3 H 2 


6 o 6 ^ 5i5 


5.8 

1.6 

0.95 

1.68 + 0.07 

217822.148 

C-C 3 H 2 


6 i 6 ^ 5o5 






217880.232 

U 





1.4 

0.44 

0.66 + 0.02 

217940.046 

C-C 3 H 2 


5i4 ^ 423 


5.7 

1.5 

0.51 

0.83 + 0.02 

218160.456 

C-C 3 H 2 


524 ^ 4i3 


5.7 

1.5 

0.23 

0.36 + 0.02 

218222.192 

H 2 CO 


3o3 ^ 2o2 


5.9 

2.4 

6.27 

16.09 + 0.03 

218324.723 

HC 3 N 


7 = 24 ^ 23 


6.1 

1.4 

0.55 

0.81+0.02 

218440.063 

CH 3 OH 


42 ^ 3i,E 


5.8 

2.1 

0.88 

1.97 + 0.03 

218475.632 

H 2 CO 


322 ^ 221 


5.9 

2.2 

1.52 

3.61 + 0.03 

218732.732 

C-C 3 H 2 


0 

t 


5.8 

1.6 

0.05 

0.11 + 0.02 

218732.732 

C-C 3 H 2 


726 ^ 7i7 






218760.066 

H 2 CO 


321 ^ 220 


5.9 

2.2 

1.35 

3.22 + 0.03 

219560.354 



J = 2^l 





20.34 + 0.04 

219798.274 

HNCO 


Ol 

T 

0 

0 

0 


5.8 

2.3 

0.13 

0.29 + 0.02 

219908.525 

WCO 


3i2 ^ 2 ii 


5.9 

2.1 

0.24 

0.54 + 0.02 
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219934.040 


A = 2 ^ 1,7 = 5/2 ^ 3/2, F = 2 ^ 1 


4.5 

0.03 

0.14 + 0.02 

219934.820 


A = 2 ^ 1,7 = 5/2 ^ 3/2, F = 3^2 





219949.442 

SO 

65 ^ 54 

5.7 

2.6 

1.96 

5.66 + 0.03 

220078.561 

CH 30 H 

80 ^ 7i, E 

6.6 

3.5 

0.06 

0.19 + 0.02 

220177.569 

H 2 CCO 

111,11 ^ 10 i,io 

5.8 

0.8 

0.10 

0.12 + 0.02 

220188.265 

u 



0.5 

0.10 

0.04 + 0.01 

220398.684 


7 = 2^1 




59.75 + 0.05 

220743.011 

CH 3 CN 

12 i ^ 111 

6.1 

1.5 

0.07 

0.12 + 0.02 

220747.261 

CH 3 CN 

12 o ^ llo 

6.1 

2.8 

0.05 

0.11 + 0.02 

221191.665 

D 2 CO 

4i4 ^ 3i3 

6.0 

1.5 

0.13 

0.20 + 0.02 

221965.220 

S 02 

111,11 ^ 10 o,io 

5.3 

2.5 

0.23 

0.60 + 0.02 

222128.815 

CH 3 CCH 

133 ^ 123 

6.0 

1.2 

0.11 

0.14 + 0.02 

222150.010 

CH 3 CCH 

132 122 

5.9 

1.3 

0.19 

0.25 + 0.01 

222162.730 

CH 3 CCH 

13i —> 12 i 

5.9 

1.5 

0.31 

0.47 + 0.02 

222166.971 

CH 3 CCH 

13o ^ 12 o 

5.9 

1.6 

0.36 

0.58 + 0.02 

223756.396 

u 



0.6 

0.08 

0.06 + 0.02 

224327.250 

H 2 CCO 

111,10 ^ 10 i ,9 

5.8 

1.5 

0.09 

0.14 + 0.02 

224714.187 


7 = 2^1 

5.2 

2.2 

2.35 

6.00 + 0.03 

225153.704 

S 02 

132,12 ^ 13i,13 

5.2 

2.3 

0.06 

0.15 + 0.02 

225159.616 

u 



1.5 

0.05 

0.07 + 0.02 

225697.775 

H 2 CO 

3i2 ^ 2 ii 

5.8 

2.5 

6.59 

19.34 + 0.04 

226004.854 

u 



3.5 

0.09 

0.31 + 0.03 

226287.418 

CN 

A = 2^ l,7 = 3/2^3/2,F= 1/2^ 1/2 

5.6 

2.1 

0.15 

0.32 + 0.02 

226298.943 

CN 

A = 2 ^ 1,7 = 3/2 ^ 3/2, F = 1/2 ^ 3/2 

5.2 

2.8 

0.11 

0.29 + 0.02 

226303.037 

CN 

A = 2^ 1,7 = 3/2 ^ 3/2,F = 3/2 ^ 1/2 

5.6 

2.0 

0.14 

0.29 + 0.02 

226314.540 

CN 

A = 2 ^ 1,7 = 3/2 ^ 3/2, F = 3/2 ^ 3/2 

5.6 

1.8 

0.33 

0.63 + 0.03 

226332.499 

CN 

A = 2 ^ 1,7 = 3/2 ^ 3/2, F = 3/2 ^ 5/2 

5.6 

2.0 

0.14 

0.31+0.02 

226341.930 

CN 

A = 2 ^ 1,7 = 3/2 ^ 3/2, F = 5/2 ^ 3/2 

5.4 

1.6 

0.19 

0.32 + 0.02 

226359.871 

CN 

A = 2 ^ 1,7 = 3/2 ^ 3/2, F = 5/2 ^ 5/2 

5.6 

1.9 

0.77 

1.66 + 0.03 

226551.622 

CH 3 CHO 

12 o,i 2 ^ llo,ii. E 

5.2 

2.6 

0.06 

0.17 + 0.03 

226592.725 

CH 3 CHO 

12 o,i 2 ^ llo,ii. A 

5.6 

2.0 

0.04 

0.08 + 0.02 

226616.571 

CN 

A = 2 ^ 1,7 = 3/2 ^ 1/2, F = 1/2 ^ 3/2 

5.6 

2.1 

0.16 

0.37 + 0.03 

226632.190 

CN 

A = 2^ 1,7 = 3/2^ 1/2,F = 3/2^3/2 

5.7 

2.1 

1.06 

2.40 + 0.03 

226659.558 

CN 

A = 2 ^ 1,7 = 3/2 ^ 1/2, F = 5/2 ^ 3/2 

5.8 

2.3 

2.44 

6.53 + 0.04 

226663.693 

CN 

A = 2^ 1,7 = 3/2^ 1/2,F= 1/2^ 1/2 

5.8 

2.3 

0.96 

2.33 + 0.03 

226679.311 

CN 

A = 2 ^ 1,7 = 3/2 ^ 1/2, F = 3/2 ^ 1/2 

5.7 

2.1 

1.22 

2.80 + 0.03 

226823.149 

U 



2.7 

0.06 

0.15 + 0.02 

226838.423 

u 



2.7 

0.03 

0.08 + 0.02 

226843.364 

u 



2.3 

0.07 

0.17 + 0.03 

226874.191 

CN 

A = 2 ^ 1,7 = 5/2 ^ 3/2, F = 5/2 ^ 3/2 


3.9 

3.50 

15.14 + 0.05 

226874.781 

CN 

A = 2 ^ 1,7 = 5/2 ^ 3/2, F = 7/2 ^ 5/2 





226875.896 

CN 

A = 2 ^ 1,7 = 5/2 ^ 3/2, F = 3/2 ^ 1/2 





226887.420 

CN 

A = 2 ^ 1,7 = 5/2 ^ 3/2, F = 3/2 ^ 3/2 

5.7 

2.0 

0.72 

1.55 + 0.03 

226892.128 

CN 

A = 2 ^ 1,7 = 5/2 ^ 3/2, F = 5/2 ^ 5/2 

5.6 

2.0 

0.84 

1.80 + 0.03 

227169.138 

C-C 3 H 2 

432 ^ 321 

5.6 

1.5 

0.50 

0.80 + 0.03 

227418.905 

HC 3 N 

7 = 25 ^ 24 

6.0 

1.4 

0.34 

0.47 + 0.02 

228910.481 

DNC 

7 = 3^2 

5.7 

1.5 

0.67 

1.02 + 0.02 

229758.756 

CH 3 OH 

8 _i ^ 7o, E 

5.3 

2.2 

0.23 

0.54 + 0.02 

229775.966 

U 



1.1 

0.05 

0.06 + 0.01 

230538.000 

CO 

7 = 2^1 




381.24 + 0.06 

231220.685 

I^CS 

7 = 5^4 

5.2 

2.0 

0.23 

0.48 + 0.02 

231410.234 

D 2 CO 

4o4 ^ 3o3 

5.5 

1.8 

0.14 

0.24 + 0.02 

233455.410 

u 



1.9 

0.06 

0.11 + 0.03 

235151.720 

S 02 

422 ^ 3i3 

5.0 

2.6 

0.08 

0.17 + 0.02 

235774.099 

CH 3 CHO 

33,1 ^ 32,2, A 

5.3 

1.0 

0.05 

0.05 + 0.01 

235816.499 

u 



0.9 

0.06 

0.05 + 0.02 

236102.086 

D 2 CO 

422 ^ 321 

5.8 

1.7 

0.10 

0.16 + 0.02 

236512.789 

HC 3 N 

7 = 26 ^ 25 

5.9 

2.0 

0.11 

0.22 + 0.02 

236727.020 

H 2 CS 

111 6 i 6 

5.6 

1.9 

0.16 

0.34 + 0.02 

237068.833 

S 02 

123,9 ^ 122,10 

5.8 

4.8 

0.04 

0.10 + 0.02 

237433.168 

u 



3.4 

0.09 

0.32 + 0.03 

237998.265 

u 



2.7 

0.10 

0.26 + 0.03 

238609.987 

u 



2.0 

0.16 

0.32 + 0.02 

239211.215 

CH 3 CCH 

143 ^ 133 

5.9 

1.2 

0.11 

0.14 + 0.02 

239234.034 

CH 3 CCH 

142 ^ 132 

5.8 

1.5 

0.13 

0.20 + 0.02 

239247.728 

CH 3 CCH 

14i ^ 13i 

5.8 

1.5 

0.27 

0.40 + 0.02 

239252.294 

CH 3 CCH 

14o ^ 13o 

5.9 

1.4 

0.32 

0.46 + 0.02 

239746.219 

CH 30 H 

5i ^4i, A+ 

5.6 

2.4 

0.28 

0.70 + 0.03 

240020.626 

u 



2.4 

0.07 

0.16 + 0.02 

240185.794 

H 2 CCO 

12 i,i 2 ^ 111,11 

5.9 

1.6 

0.09 

0.15 + 0.02 
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240266.872 

H 2 CS 


101 ^ 6 o 6 


5.6 

2.1 

0.13 

0.29 + 0.02 

241016.089 



7 = 5^4 


5.7 

1.8 

0.86 

1.68 + 0.03 

241615.797 

SO 2 


524 ^ 4i3 


5.5 

2.2 

0.16 

0.39 + 0.03 

241700.159 

CH 3 OH 


5o ^ 4o, E 


5.7 

2.2 

0.50 

1.12 + 0.03 

241767.234 

CH 3 OH 


5_i ^4_i,E 


5.6 

2.2 

1.27 

3.17 + 0.04 

241774.032 

HNCO 


llo,ii ^ 10 o,io 


5.7 

2.0 

0.09 

0.23 + 0.03 

241791.352 

CH 3 OH 


5o ^ 4o, A+ 


5.6 

2.2 

1.50 

3.80 + 0.05 

241832.718 

CH 3 OH 


53 ^ 43 , A+ 



1.7 

0.06 

0.12 + 0.03 

241833.106 

CH 3 OH 


5_3^4_3,A- 






241842.284 

CH 3 OH 


5-2^4_2,A- 



2.8 

0.05 

0.12 + 0.02 

241843.604 

CH 3 OH 


53 ^ 43 , E 






241879.025 

CH 3 OH 


5i ^4i,E 


5.6 

2.1 

0.25 

0.53 + 0.03 

241887.674 

CH 3 OH 


52 ^ 42, A+ 


5.6 

2.1 

0.05 

0.10 + 0.03 

241904.147 

CH 3 OH 


5-2^4_2,E 



2.4 

0.42 

1.04 + 0.03 

241904.643 

CH 3 OH 


52 ^ 42, E 






242913.610 



7 = 5^4 


5.6 

1.8 

0.19 

0.37 + 0.02 

243915.788 

CH 30 H 


5_i ^4_i,A- 


5.5 

2.2 

0.31 

0.71 + 0.03 

244048.504 

H 2 CS 


7i6 ^ 615 


5.7 

2.0 

0.20 

0.43 + 0.02 

244222.150 

C-C 3 H 2 


321 ^ 2 i 2 


5.4 

1.5 

0.15 

0.22 + 0.02 

244254.218 

S 02 


14o,i4 ^ 13i,i3 


5.2 

1.6 

0.13 

0.21 + 0.02 

244712.269 

H 2 CCO 


12 i,ii ^ 111,10 


5.8 

2.1 

0.06 

0.14 + 0.02 

244935.557 

cs 


7 = 5^4 


5.8 

2.3 

8.75 

22.76 + 0.04 

245532.752 

D 2 CO 


4i3 ^ 3i2 


5.8 

1.8 

0.06 

0.11 + 0.02 

337061.050 



7 = 3^2 


5.5 

1.7 

2.64 

5.19 + 0.04 

337396.459 



7 = 7^6 


5.8 

1.8 

0.18 

0.33 + 0.02 

338072.030 

u 





2.2 

0.06 

0.13 + 0.02 

338083.195 

H 2 CS 


10 i,io ^ 9 i ,9 


5.9 

1.7 

0.09 

0.16 + 0.02 

338124.488 

CH 30 H 


7o ^ 6 o» E 


5.7 

1.7 

0.15 

0.28 + 0.02 

338344.588 

CH 30 H 


7_i ^ 6 _i,E 


5.7 

1.8 

0.43 

0.83 + 0.03 

338408.698 

CH 30 H 


7o ^ 6 o» A+ 


5.7 

1.8 

0.48 

0.92 + 0.03 

338614.936 

CH 30 H 


7i ^ 61 , E 


5.7 

2.3 

0.09 

0.22 + 0.03 

338721.693 

CH 30 H 


72 ^ 62 , E 


5.7 

1.9 

0.17 

0.35 + 0.02 

339341.459 

so 


33 ^ 32 


5.8 

1.5 

0.10 

0.18 + 0.02 

339475.904 

CN 

N = 

3^2,7 = 5/2^5/2,F = 

5/2 ^ 5/2 

5.7 

1.7 

0.09 

0.17 + 0.02 

339516.635 

CN 


3 ^2,7 = 5/2 ^5/2, F = 

7/2 ^ 7/2 

5.7 

1.5 

0.15 

0.24 + 0.02 

340008.126 

CN 

N = 

3^2,7 = 5/2^3/2,F = 

5/2 ^ 5/2 

5.7 

1.7 

0.20 

0.34 + 0.02 

340019.626 

CN 


3^2,7 = 5/2^3/2,F = 

3/2 ^ 3/2 

5.8 

1.9 

0.19 

0.36 + 0.02 

340031.549 

CN 

N = 

3^2,7 = 5/2^3/2,F = 

7/2 ^ 5/2 

5.7 

2.2 

1.19 

5.33 + 0.04 

340035.408 

CN 


3 ^ 2,7 = 5/2^ 3/2, F = 

5/2 ^ 3/2 





340035.408 

CN 

N = 

3^2,7 = 5/2^3/2,F = 

3/2 ^ 1/2 





340247.770 

CN 


3 ^2,7 = 7/2 ^5/2, F = 

7/2 ^ 5/2 

5.7 

2.6 

2.52 

7.36 + 0.04 

340247.770 

CN 

N = 

3^2,7 = 7/2^5/2,F = 

9/2 ^ 7/2 





340248.544 

CN 


3 ^ 2,7 = 7/2^ 5/2, F = 

5/2 ^ 3/2 





340261.773 

CN 

N = 

3^2,7 = 7/2^5/2,F = 

5/2 ^ 5/2 

5.8 

1.8 

0.19 

0.33 + 0.02 

340264.949 

CN 


3 ^2,7 = 7/2 ^5/2, F = 

7/2 ^ 7/2 

5.6 

1.4 

0.22 

0.33 + 0.02 

340630.692 

HCi^O+ 


7 = 4^3 


5.8 

1.8 

0.13 

0.24 + 0.02 

340714.155 

SO 


78 ^ 67 


5.8 

1.9 

0.60 

1.34 + 0.03 

362045.753 

DCN 


7 = 5^4 


6.0 

1.3 

0.27 

0.34 + 0.05 

362630.303 

HNC 


7 = 4^3 


5.9 

1.7 

5.18 

9.63 + 0.09 

362736.048 

H 2 CO 


5o5 ^ 4o4 


6.0 

1.8 

3.47 

6.77 + 0.09 

363739.868 

CH 3 OH 


72 ^ 6 i,E 


5.7 

2.0 

0.19 

0.43 + 0.07 

363945.894 

H 2 CO 


524 ^ 423 


5.9 

1.8 

0.75 

1.37 + 0.07 

364275.141 

H 2 CO 


533 ^ 432 


5.9 

2.0 

0.31 

0.70 + 0.07 

364288.884 

H 2 CO 


532 ^ 431 


5.7 

1.8 

0.29 

0.51+0.07 

365363.428 

H 2 CO 


523 ^ 422 


5.8 

2.4 

0.49 

1.33 + 0.10 


Notes. Rest frequency. For the unidentified lines, this is calculated assuming a source velocity of 5.7 km s 


Table B.2. Observed line strengths towards CrA-46. 


Frequency Molecule Transition vlsr Av Tnib f T^\ydv 

[MHz] [kms-i] [kms-i] [K] [Kkms-^] 


No detections 
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Table B.3. Observed line strengths towards CrA-3. 


Frequency 

[MHz] 

Molecule 

Transition 

Vlsr 
[ km s"^] 

Av 

[km s"^] 

Tmh 

[K] 

/ Tmhdv 
[Kkms-'l 

219560.354 

c'*o 

7 = 2^1 

6.0 

1.1 

0.63 

0.69 ± 0.03 

337061.050 

c‘^o 

7 = 3^2 

5.8 

0.3 

0.26 

0.06 + 0.03 


Table B.4. Observed line strengths towards CrA-5. 


Frequency Molecule Transition vlsr Av Tmb f Tvcib^V 

[MHz] [kms-i] [kms-i] [K] [Kkms-^] 

219560.354 C^^O / = 2 ^ 1 6.2 0.7 1.96 1.40 + 0.02 


Table B.5. Observed line strengths towards LS-RCrAl. 


Frequency Molecule Transition Vlsr Av ^mb f 

[MHz] [kms-^] [kms-i] [K] [Kkms"^] 

219560.354 C^^O 7 = 2^ 1 6.0 1.4 0.86 1.22 + 0.04 


Table B.6. Observed line strengths towards Haas4. 


Frequency 

Molecule 

Transition 

Vlsr 

Av 

7’mb 

/ T^dv 

[MHz] 



[km s-‘] 

[km s"^] 

[K] 

[Kkms-‘] 

216112.582 

DCO+ 

7 = 3^2 

5.6 

0.9 

0.61 

0.55 + 0.02 

216278.756 

C-C 3 H 2 

3o3 ^ 221 

5.4 

0.5 

0.24 

0 . 11 + 0.02 

217104.919 

SiO 

7 = 5^4 

10.4 

5.6 

0.15 

1.08 + 0.06 

217822.148 

C-C 3 H 2 

6o6 ^ 5i5 

5.4 

0.6 

0.12 

0.07 + 0.02 

217822.148 

C-C 3 H 2 

6i6 ^ 5o5 





217940.046 

C-C 3 H 2 

5i4 ^ 423 

5.7 

0.9 

0.11 

0.12 + 0.03 

218222.192 

H 2 CO 

3o3 ^ 2o2 

5.5 

1.0 

1.18 

1.84 + 0.05 

219560.354 


7 = 2^1 

5.5 

1.3 

9.79 

14.62 + 0.04 

219949.442 

SO 

65 ^ 54 

5.5 

0.9 

1.15 

2.30 + 0.06 

337061.050 


7 = 3^2 

5.4 

1.3 

0.84 

1.20 + 0.04 

340247.770 

CN 

A = 3 ^ 2,7 = 7/2 ^ 5/2, F = 7/2 ^ 5/2 

5.4 

0.6 

0.20 

0.12 + 0.02 

340247.770 

CN 

A = 3 ^ 2,7 = 7/2 ^ 5/2, F = 9/2 ^ 7/2 






Table B.7. Observed line strengths towards IRS2. 


Frequency Molecule Transition vlsr Av 7’mb f T^dv 

[MHz] [kms-‘] [kms-'l [K] [Kkms-‘] 


216112.582 

DCO+ 


7 = 3^2 


6.0 

0.4 

0.87 

0.37 + 0.01 

218222.192 

H 2 CO 


3o3 ^ 2o2 


6.0 

0.6 

1.97 

1.29 + 0.01 

218440.063 

CH 3 OH 


42^3i,E 


6.0 

0.4 

0.16 

0.06 + 0.01 

218475.632 

H 2 CO 


322 ^ 221 


6.0 

0.6 

0.07 

0.05 + 0.01 

218760.066 

H 2 CO 


321 ^ 220 


6.2 

0.9 

0.07 

0.07 + 0.01 

219560.354 



7 = 2^1 


6.2 

1.0 

7.72 

8.47 + 0.02 

219908.525 

H^^CO 


3i 2 ^ 2ii 


6.0 

0.5 

0.07 

0.04 + 0.01 

219949.442 

SO 


65 ^ 54 


6.0 

0.4 

2.47 

1.22 + 0.01 

337061.050 



7 = 3^2 


6.1 

1.0 

0.97 

1 . 11 + 0.02 

339341.459 

SO 


33^32 


6.1 

0.4 

0.19 

0.08 + 0.01 

340247.770 

CN 

A = 

3 ^2,7 = 7/2^5/2,F 

= 7/2 ^ 5/2 

6.1 

1.9 

0.09 

0.15 + 0.02 

340247.770 

CN 

A = 

3 ^ 2,7 = 7/2^ 5/2, F 

= 9/2 ^ 7/2 
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Table B.8. Observed line strengths towards IRS5A. 


Frequency 

[MHz] 

Molecule 


Transition 


VlSR 
[km s"^] 

Av 

[km s"^] 

^mb 

[K] 

f Tn,bdv 
[Kkms-*] 

216112.582 

DCO+ 


7 = 3^2 


5.7 

0.7 

0.99 

0.81+0.02 

216278.756 

C-C3H2 


3o3 ^ 221 


5.6 

0.9 

0.38 

0.39 + 0.02 

216373.320 

C2D 

A = 3 ^ 2,/ = 

7/2^5/2,F = 

5/2 ^ 3/2 

5.9 

1.5 

0.16 

0.26 + 0.02 

216373.320 

C2D 

A = 3 ^ 2,/ = 

7/2^ 5/2, F = 

7/2 ^ 5/2 





216428.320 

C2D 

A = 3 ^ 2,/ = 

5/2^3/2,F = 

7/2 ^ 5/2 

5.5 

1.0 

0.12 

0.12 + 0.01 

216428.320 

C2D 

A = 3 ^ 2,/ = 

512 ^ 312, 

5/2 ^ 3/2 





217238.538 

DCN 


7 = 3^2 


5.6 

1.0 

0.18 

0.20 + 0.02 

217467.150 

'^CN 

N = 1,/ = 5/2 

^ 3/2, Fi = 3 - 

-^2,F = 4^3 

5.4 

0.8 

0.11 

0.09 + 0.01 

217467.150 

‘^CN 

N = 1,/ = 5/2 

^ 3/2,Fi = 3 - 

2, F = 3 ^ 2 





217822.148 

C-C3H2 


6o6 ^ 5i5 


5.6 

1.0 

0.34 

0.37 + 0.02 

217822.148 

C-C3H2 


6i6 ^ 5o5 






217940.046 

C-C3H2 


5i4 ^ 423 


5.6 

1.0 

0.20 

0.21 + 0.02 

218222.192 

H2CO 


%3 2o2 


5.6 

1.3 

3.26 

4.87 + 0.03 

218440.063 

CH3OH 


42^3i,E 


5.6 

1.1 

0.20 

0.26 + 0.02 

218475.632 

H2CO 


322 ^ 221 


5.5 

1.4 

0.37 

0.58 + 0.02 

218760.066 

H2CO 


321 ^ 220 


5.5 

1.4 

0.38 

0.54 + 0.02 

219560.354 

c‘*o 


7 = 2^1 


5.7 

1.1 

12.72 

16.13 + 0.03 

219908.525 

H'^CO 


3i2 ^ 2ii 


5.6 

1.3 

0.09 

0.13 + 0.02 

219949.442 

SO 


65 ^ 54 


5.7 

0.7 

6.05 

4.80 + 0.02 

337061.050 



7 = 3^2 


5.6 

1.2 

2.03 

2.65 + 0.03 

338344.588 

CH30H 

7-1 ^ 6 _„E 


5.5 

1.6 

0.10 

0.16 + 0.03 

338408.698 

CH30H 

2o ^ 60 , A+ 


5.4 

1.6 

0.14 

0.23 + 0.03 

339341.459 

so 


33 ^ 32 


5.7 

0.5 

0.73 

0.40 + 0.02 

340008.126 

CN 

A = 3 ^ 2,/ = 

5/2^3/2,F = 

5/2 ^ 5/2 

5.4 

1.5 

0.14 

0.23 + 0.03 

340019.626 

CN 

A = 3 ^2,/ = 

5/2^3/2,F = 

3/2 ^ 3/2 

5.6 

1.1 

0.15 

0.21 + 0.03 

340031.549 

CN 

N = 3 ^ 2, J = 

5/2^3/2,F = 

7/2 ^ 5/2 

5.4 

1.7 

0.69 

1.17 + 0.03 

340035.408 

CN 

A = 3 ^2,/ = 

5/2^3/2,F = 

3/2 ^ 1/2 

5.4 

1.6 

0.69 

1.11+0.03 

340035.408 

CN 

N = 3^2,J = 

5/2 ^ 3/2, F = 

5/2 ^ 3/2 





340247.770 

CN 

A = 3 ^2,/ = 

7/2^5/2,F = 

7/2 ^ 5/2 

5.3 

2.2 

1.42 

3.36 + 0.04 

340247.770 

CN 

N = 3 ^ 2, J = 

7/2^ 5/2, F = 

9/2 ^ 7/2 





340261.773 

CN 

A = 3 ^2,7 = 

7/2^ 5/2, F = 

5/2 ^ 5/2 

5.7 

2.5 

0.10 

0.16 + 0.02 

340264.949 

CN 

A = 3 ^2,7 = 

7/2^ 5/2, F = 

7/2 ^ 7/2 

5.5 

1.2 

0.20 

0.23 + 0.03 

340630.692 

HCi^O+ 


7 = 4^3 


5.5 

0.6 

0.13 

0.07 + 0.02 

340714.155 

SO 


78 ^ 67 


5.5 

0.5 

0.52 

0.27 + 0.02 

362630.303 

HNC 


7 = 4^3 


5.6 

1.5 

1.59 

2.41 + 0.06 

362736.048 

H 2 CO 


5o5 ^ 4o4 


5.6 

1.2 

0.74 

0.93 + 0.05 
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Table B.9. Observed line strengths towards IRS5N. 


Frequency 

Molecule 

Transition 


VlSR 

Av 

Tmh 

f Tn,bdv 

[MHz] 




[km s"^] 

[km s"^] 

[K] 

[Kkms-*] 

216112.582 

DCO+ 

7 = 3^2 


5.7 

0.9 

0.80 

0.78 + 0.02 

216278.756 

C-C3H2 

3o3 ^ 221 


5.7 

1.3 

0.26 

0.35 + 0.02 

216373.320 

C2D 

A = 3 ^2,/ = 7/2^5/2,F = 

5/2 ^ 3/2 

5.9 

1.6 

0.17 

0.29 + 0.02 

216373.320 

C2D 

A = 3 ^2,/ = 7/2^5/2,F = 

7/2 ^ 5/2 





216428.320 

C2D 

A = 3 ^2,/ = 5/2^3/2,F = 

7/2 ^ 5/2 

5.5 

1.6 

0.10 

0.15 + 0.02 

216428.320 

C2D 

A = 3 ^2,/ = 5/2^3/2,F = 

5/2 ^ 3/2 





217238.538 

DCN 

7 = 3^2 


5.8 

1.5 

0.17 

0.26 + 0.02 

217467.150 

'^CN 

1,7 = 5/2^ 3/2,Fi = 3- 

-^2,F = 4^3 

5.6 

1.6 

0.08 

0.12 + 0.02 

217467.150 

‘^CN 

N = 1,7 = 5/2^ 3/2, Fi = 3^2,F = 3^2 





217822.148 

C-C3H2 

6 o 6 ^ 5i5 


5.6 

1.2 

0.27 

0.36 + 0.02 

217822.148 

C-C3H2 

6 i 6 ^ 5o5 






217940.046 

C-C3H2 

5i4 ^ 423 


5.7 

1.4 

0.16 

0.23 + 0.02 

218160.456 

C-C3H2 

524 ^ 4i3 


5.7 

1.1 

0.07 

0.09 + 0.02 

218222.192 

H2CO 

3o3 ^ 2o2 


5.7 

1.4 

3.26 

4.95 + 0.03 

218440.063 

CH3OH 

42^3i,E 


5.7 

1.2 

0.23 

0.30 + 0.02 

218475.632 

H2CO 

322 ^ 221 


5.6 

1.4 

0.35 

0.50 + 0.02 

218760.066 

H2CO 

321 ^ 220 


5.6 

1.3 

0.37 

0.54 + 0.02 

219560.354 

c'*o 

7 = 2^1 


5.7 

1.1 

12.39 

15.71 + 0.03 

219908.525 

HfCO 

3i2 ^ 2ii 


5.8 

1.1 

0.12 

0.14 + 0.02 

219934.040 

C'^N 

A = 2 ^ 1,7 = 5/2 ^ 3/2,F 

= 2^1 

4.3 

2.7 

0.04 

0.11+0.03 

219934.820 

C'^N 

N = 2^ l,7 = 5/2^3/2,F 

= 3^2 





219949.442 

SO 

65 ^ 54 


5.7 

0.7 

3.77 

3.00 + 0.02 

337061.050 


7 = 3^2 


5.7 

1.2 

1.94 

2.53 + 0.03 

338344.588 

CH30H 

7-1 ^6.1, E 


5.5 

1.2 

0.15 

0.18 + 0.02 

338408.698 

CH30H 

2o ^ 60, A+ 


5.6 

1.4 

0.16 

0.25 + 0.03 

339341.459 

so 

33 ^ 32 


5.7 

0.5 

0.29 

0.15 + 0.02 

339516.635 

CN 

7V = 3 ^2,7 = 5/2^5/2,F = 

7/2 ^ 7/2 

5.4 

1.1 

0.10 

0.14 + 0.02 

340008.126 

CN 

A = 3 ^2,7 = 5/2^3/2,F = 

5/2 ^ 5/2 

5.4 

0.9 

0.15 

0.18 + 0.02 

340019.626 

CN 

A = 3 ^2,7 = 5/2^3/2,F = 

3/2 ^ 3/2 

5.6 

1.5 

0.13 

0.20 + 0.02 

340031.549 

CN 

A = 3 ^2,7 = 5/2^3/2,F = 

7/2 ^ 5/2 

5.6 

1.8 

0.72 

1.15 + 0.02 

340035.408 

CN 

A = 3 ^2,7 = 5/2^3/2,F = 

3/2 ^ 1/2 

5.6 

1.5 

0.69 

1.03 + 0.02 

340035.408 

CN 

A = 3 ^2,7 = 5/2^3/2,F = 

5/2 ^ 3/2 





340247.770 

CN 

A = 3 ^2,7 = 7/2^5/2,F = 

7/2 ^ 5/2 

5.4 

2.1 

1.44 

3.15 + 0.03 

340247.770 

CN 

A = 3 ^2,7 = 7/2^5/2,F = 

9/2 ^ 7/2 





340261.773 

CN 

A = 3 ^2,7 = 7/2^5/2,F = 

5/2 ^ 5/2 

5.5 

1.0 

0.13 

0.12 + 0.02 

340264.949 

CN 

A = 3 ^2,7 = 7/2^5/2,F = 

7/2 ^ 7/2 

5.6 

1.2 

0.16 

0.19 + 0.02 

340714.155 

SO 

78 ^ 67 


5.5 

0.6 

0.17 

0.12 + 0.02 

362150.763 

u 




0.7 

0.24 

0.15 + 0.04 

362630.303 

HNC 

7 = 4^3 


5.9 

1.5 

1.83 

2.75 + 0.07 

362736.048 

H2CO 

5o5 ^ 4o4 


5.7 

1.2 

0.77 

0.90 + 0.06 
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Table B.IO. Observed line strengths towards IRSL 


Frequency 

Molecule 


Transition 


VlSR 

Av 

Tmh 

f Tn,bdv 

[MHz] 





[km s"^] 

[km s"^] 

[K] 

[Kkms-*] 

216112.582 

DCO+ 


7 = 3^2 


5.8 

0.8 

0.44 

0.38 ± 0.02 

216278.756 

C-C3H2 


3o3 ^ 221 


5.8 

0.8 

0.15 

0.13 + 0.02 

216373.320 

C2D 

A = 3 ^ 2,/ = 

7/2^5/2,F = 

5/2 ^ 3/2 

6.0 

1.2 

0.10 

0.13 + 0.02 

216373.320 

C2D 

A = 3 ^ 2,/ = 

7/2^ 5/2, F = 

7/2 ^ 5/2 





216428.320 

C2D 

A = 3 ^ 2,/ = 

5/2^3/2,F = 

7/2 ^ 5/2 

5.8 

0.8 

0.09 

0.06 + 0.02 

216428.320 

C2D 

A = 3 ^ 2,/ = 

512 ^ 312, 

5/2 ^ 3/2 





217238.538 

DCN 


7 = 3^2 


5.8 

0.9 

0.17 

0.17 + 0.02 

217467.150 

'^CN 

N = 1,/ = 5/2 

^ 3/2, Fi = 3 - 

-^2,F = 4^3 

5.7 

2.5 

0.05 

0.15 + 0.03 

217467.150 

‘^CN 

N = 1,/ = 5/2 

^ 3/2,Fi = 3 - 

2, F = 3 ^ 2 





217822.148 

C-C3H2 


6 o6 ^ 5i5 


5.8 

0.8 

0.13 

0.10 + 0.02 

217822.148 

C-C3H2 


6 i6 ^ 5o5 






217940.046 

C-C3H2 


5i4 ^ 423 


5.6 

0.6 

0.09 

0.05 + 0.02 

218222.192 

H2CO 


%3 2 o2 


5.8 

1.2 

3.05 

4.08 + 0.03 

218440.063 

CH3OH 


42^3i,E 


5.9 

1.2 

0.30 

0.39 + 0.02 

218475.632 

H2CO 


322 ^ 221 


5.8 

1.4 

0.33 

0.54 + 0.02 

218760.066 

H2CO 


321 ^ 220 


5.8 

1.3 

0.35 

0.39 + 0.02 

219560.354 

c‘*o 


7 = 2^1 


5.8 

1.1 

11.61 

13.63 + 0.03 

219908.525 

H'^CO 


3i2 ^ 2 ii 


5.9 

0.9 

0.09 

0.09 + 0.02 

219949.442 

SO 


65 ^ 54 


5.8 

1.0 

2.20 

2.47 + 0.03 

337061.050 



7 = 3^2 


5.8 

1.1 

1.70 

1.97 + 0.03 

338344.588 

CH30H 

7-1 ^6_„E 


5.8 

1.1 

0.24 

0.26 + 0.03 

338408.698 

CH30H 

2o ^ 60, A+ 


5.8 

1.4 

0.24 

0.36 + 0.03 

340031.549 

CN 

A = 3 ^ 2,/ = 

5/2 ^ 3/2, F = 

7/2 ^ 5/2 

5.9 

1.0 

0.49 

0.52 + 0.03 

340035.408 

CN 

N = 3 ^ 2, J = 

5/2^3/2,F = 

3/2 ^ 1/2 

5.8 

1.0 

0.47 

0.48 + 0.03 

340035.408 

CN 

A = 3 ^2,/ = 

5/2^3/2,F = 

5/2 ^ 3/2 





340247.770 

CN 

N = 3 ^ 2, J = 

7/2^ 5/2, F = 

7/2 ^ 5/2 

5.7 

1.6 

0.87 

1.48 + 0.04 

340247.770 

CN 

N = 3 ^2, J = 

7/2^5/2,F = 

9/2 ^ 7/2 





340261.773 

CN 

N = 3^2,J = 

7/2^ 5/2, F = 

5/2 ^ 5/2 

5.8 

0.8 

0.11 

0.09 + 0.02 

340264.949 

CN 

A = 3 ^2,/ = 

7/2^5/2,F = 

7/2 ^ 7/2 

5.9 

0.8 

0.13 

0.10 + 0.02 

340714.155 

SO 


78 ^ 67 


5.6 

1.7 

0.32 

0.60 + 0.04 

362630.303 

HNC 


7 = 4^3 


6.0 

1.4 

0.73 

1.02 + 0.06 

362736.048 

H2CO 


5o5 ^ 4o4 


6.0 

1.7 

0.48 

0.78 + 0.06 
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Table B.ll. Observed line strengths towards IRS7A. 


Frequency 

Molecule 



Transition 



VlSR 

Av 

7’mb 

f Trnbdv 

[MHz] 







[km s"^] 

[km s"^] 

[K] 

[Kkms-‘] 

337061.050 




7 = 3^2 



5.1 

2.0 

2.39 

5.63 ± 0.04 

337396.459 




J^l ^6 



5.0 

2.4 

0.21 

0.54 ± 0.03 

337580.147 




88^7v 



4.6 

3.8 

0.07 

0.34 + 0.04 

338083.195 

H 2 CS 



10 i,io ^ 9 i ,9 



5.2 

1.4 

0.08 

0 . 11 + 0.02 

338124.488 

CH 30 H 



7o ^ 60 , E 



5.0 

2.2 

0.12 

0.24 + 0.02 

338344.588 

CH 30 H 



7-1 ^ 6 _i,E 



5.1 

2.4 

0.36 

0.87 ± 0.03 

338408.698 

CH 30 H 



2o ^ 60 , A+ 



5.1 

2.3 

0.45 

1.09 + 0.03 

338614.936 

CH 30 H 



7i ^ 6 i,E 



5.2 

2.2 

0.10 

0.19 + 0.03 

338721.693 

CH 30 H 



li 62 , E 



4.9 

2.8 

0.17 

0.40 + 0.03 

339341.459 

so 



33 ^ 32 



4.9 

2.6 

0.19 

0.52 + 0.03 

339516.635 

CN 

A = 3 


2,/ = 5/2^5/2,F = 7/2 


111 

5.1 

1.3 

0.12 

0.15 + 0.02 

339857.269 




98 ^ 87 



4.6 

2.8 

0.14 

0.42 + 0.03 

340008.126 

CN 

A = 3 


2,/ = 5/2^3/2,F = 5/2 


5/2 

5.1 

2.3 

0.16 

0.39 + 0.03 

340019.626 

CN 

A = 3 


2,/ = 5/2^3/2,F = 3/2 


3/2 

5.1 

2.0 

0.17 

0.39 + 0.03 

340031.549 

CN 

A = 3 


2,/ = 5/2^3/2,F = 7/2 


5/2 

5.1 

2.4 

1.13 

5.67 + 0.05 

340035.408 

CN 

A = 3 


2,/ = 5/2^3/2,F = 5/2 


3/2 





340035.408 

CN 

A = 3 


2,/ = 5/2^3/2,F = 3/2 


1/2 





340247.770 

CN 

A = 3 


2,/ = 7/2^5/2,F = 7/2 


5/2 

4.9 

2.8 

2.64 

8.14 + 0.04 

340247.770 

CN 

A = 3 


2,/ = 7/2^5/2,F = 9/2 


7/2 





340261.773 

CN 

A = 3 


2,/ = 7/2^5/2,F = 5/2 


5/2 

5.3 

1.9 

0.14 

0.31 + 0.02 

340264.949 

CN 

A = 3 


2,/ = 7/2^5/2,F = 7/2 


7/2 

5.0 

2.0 

0.16 

0.36 + 0.03 

340630.692 

HCi^O+ 



7 = 4^3 



5.1 

1.3 

0.19 

0.29 + 0.03 

340714.155 

SO 



78 ^ 67 



4.7 

3.0 

1.18 

3.66 + 0.04 

362045.753 

DCN 



7 = 5^4 



5.8 

2.6 

0.19 

0.53 + 0.08 

362630.303 

HNC 



7 = 4^3 



5.5 

2.3 

3.13 

7.57 + 0.09 

362736.048 

H 2 CO 



5o5 ^ 4o4 



5.4 

2.5 

2.30 

6.39 + 0.10 

363945.894 

H 2 CO 



524 ^ 423 



5.4 

2.4 

0.54 

1.35 + 0.08 

364275.141 

H 2 CO 



533 ^ 432 



5.3 

2.5 

0.21 

0.61 + 0.10 

364288.884 

H 2 CO 



532 ^ 431 



5.6 

2.7 

0.20 

0.59 + 0.10 
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Table B.12. Observed line strengths towards CrA-24, total flux (both souree and outflow components). 


Frequency 

[MHz] 

Molecule 


Transition 


VlSR 
[km s"^] 

Av 

[km s"^] 

Tmb 

[K] 

f T„tbdv 
[Kkms-‘] 

216112.582 

DCO+ 


7 = 3^2 


4.9 

1.2 

0.26 

0.31+0.03 

216278.756 

C-C 3 H 2 


3o3 ^ 221 


4.8 

0.9 

0.21 

0.24 + 0.02 

216373.320 

C 2 D 


3 ^ 2,7 = 7/2^ 5/2, F 

= 5/2 ^ 3/2 

5.2 

1.5 

0.10 

0.19 + 0.04 

216373.320 

C 2 D 


3 ^ 2,7 = 7/2^ 5/2, F 

= 7/2 ^ 5/2 





217104.919 

SiO 


7 = 5^4 


1.9 

4.5 

0.10 

0.49 + 0.05 

217238.538 

DCN 


7 = 3^2 


4.7 

1.1 

0.24 

0.27 + 0.03 

217822.148 

C-C 3 H 2 


6 o 6 ^ 5i5 


4.8 

1.3 

0.27 

0.37 + 0.04 

217822.148 

C-C 3 H 2 


6i6 ^ 5o5 






217940.046 

C-C 3 H 2 


5i4 ^ 423 


5.0 

1.0 

0.16 

0.19 + 0.04 

218222.192 

H 2 CO 


3o3 ^ 2o2 





7.99 + 0.07 

218440.063 

CH 3 OH 


42^3i,E 


4.8 

2.2 

0.17 

0.47 + 0.06 

218475.632 

H 2 CO 


322 ^ 221 





1.47 + 0.06 

218760.066 

H 2 CO 


321 ^ 220 





1.63 + 0.06 

219560.354 

c‘*o 


7 = 2^1 


5.1 

1.5 

9.45 

17.30 + 0.07 

219908.525 

H'^CO 


3i2 ^ 2ii 


4.6 

1.6 

0.10 

0.16 + 0.04 

219949.442 

SO 


65 ^ 54 





2.40 + 0.06 

337061.050 



7 = 3^2 


4.9 

1.7 

1.60 

3.11 + 0.05 

338344.588 

CH 30 H 


7-1 ^ 6 _i,E 





0.60 + 0.04 

338408.698 

CH 30 H 


2o ^ 60, A-i- 





0.73 + 0.03 

340008.126 

CN 


3 ^2,7 = 5/2^3/2,F 

= 5/2 ^ 5/2 

4.6 

1.1 

0.11 

0.13 + 0.03 

340031.549 

CN 


3 ^ 2,7 = 5/2^ 3/2, F 

= 7/2 ^ 5/2 

4.7 

1.7 

0.43 

0.76 + 0.04 

340035.408 

CN 


3 ^2,7 = 5/2^3/2,F 

= 3/2 ^ 1/2 

4.7 

1.9 

0.40 

0.77 + 0.04 

340035.408 

CN 


3 ^ 2,7 = 5/2^ 3/2, F 

= 5/2 ^ 3/2 





340247.770 

CN 


3 ^2,7 = 7/2^5/2,F 

= 7/2 ^ 5/2 

4.6 

1.7 

1.12 

2.49 + 0.06 

340247.770 

CN 


3 ^ 2,7 = 7/2^ 5/2, F 

= 9/2 ^ 7/2 





340714.155 

SO 


7g 67 


5.0 

1.7 

0.12 

0.28 + 0.05 

362630.303 

HNC 


7 = 4^3 


4.6 

1.3 

1.98 

2.67 + 0.07 

362736.048 

H 2 CO 


5o5 ^ 4o4 





2.43 + 0.09 

363945.894 

H 2 CO 


524 ^ 423 


4.7 

1.3 

0.25 

0.32 + 0.07 


Table B.13. Observed line strengths towards CrA-24, on-source velocity component. The integrated intensities are calculated from Gaussian fits. 


Frequency 

Molecule 


Transition 


VlSR 

Av 

Tmb 

/ Tra>dv 

[MHz] 





[km s"^] 

[km s"^] 

[K] 

[Kkms-‘] 

217104.919 

SiO 


7 = 5^4 


4.6 

2.4 

0.06 

0.17 + 0.03 

218222.192 

H2CO 


%3 2 o 2 


5.0 

1.8 

2.61 

5.09 + 0.03 

218440.063 

CH3OH 


42^3i,E 


5.1 

1.5 

0.17 

0.27 + 0.03 

218475.632 

H2CO 


322 ^ 221 


5.0 

1.9 

0.42 

0.84 + 0.03 

218760.066 

H2CO 


321 ^ 220 


4.9 

1.9 

0.46 

0.92 + 0.03 

219908.525 

H'^CO 


3i2 ^ 2 ii 


4.6 

1.6 

0.10 

0.16 + 0.03 

219949.442 

SO 


65 ^ 54 


5.1 

1.7 

0.87 

1.58 + 0.03 

338344.588 

CH3OH 


7-1 ^6_1,E 


5.1 

2.0 

0.30 

0.63 + 0.03 

338408.698 

CH3OH 


7o ^ 60 , A+ 


4.9 

1.7 

0.42 

0.77 + 0.03 

340247.770 

CN 

N 

= 3 ^ 2,7 = 7/2^ 5/2, F 

= 7/2 ^ 5/2 

4.6 

1.7 

1.13 

2.00 + 0.03 

340247.770 

CN 

N 

= 3 ^2,7 = 7/2^5/2,F 

= 9/2 ^ 7/2 





340714.155 

SO 


Is 67 


5.0 

1.7 

0.12 

0.22 + 0.03 

362736.048 

H2CO 


5o5 ^ 4o4 


4.9 

1.7 

1.13 

2.07 + 0.05 

363945.894 

H2CO 


524 ^ 423 


4.7 

1.2 

0.26 

0.32 + 0.04 



Table B.14. Observed line strengths 
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towards CrA-24, outflow velocity component. The integrated intensities 

are calculated from Gaussian fits. 

Frequency 

Molecule 


Transition 


VlSR 

Av 

Tmb 

f T^hdv 

[MHz] 





[km s"^] 

[km s"^] 

[K] 

[Kkms-‘] 

217104.919 

SiO 


7 = 5^4 


1.2 

2.2 

0.13 

0.30 ± 0.03 

218222.192 

H 2 CO 


3o3 ^ 2o2 


2.0 

3.2 

0.75 

2.57 + 0.04 

218440.063 

CH 3 OH 


42^3i,E 


2.7 

3.5 

0.06 

0.23 ± 0.04 

218475.632 

H 2 CO 


322 ^ ^21 


2.2 

3.7 

0.18 

0.71 + 0.04 

218760.066 

H 2 CO 


321 ^ ^20 


1.6 

3.7 

0.18 

0.73 ± 0.04 

219908.525 

H'^CO 


3i2 ^ 2ii 


1.3 

1.1 

0.05 

0.06 ± 0.02 

219949.442 

SO 


65 ^ 54 


1.9 

2.8 

0.25 

0.74 + 0.03 

338344.588 

CH 3 OH 


7-1 ^ 6 _ 1 ,E 


2.0 

2.3 

0.13 

0.33 + 0.03 

338408.698 

CH 3 OH 


7o ^ 60 , A+ 


1.6 

2.5 

0.14 

0.38 ± 0.03 

340247.770 

CN 


3 ^ 2,7 = 7/2^ 5/2, F 

= 7/2 ^ 5/2 

1.6 

2.9 

0.11 

0.34 + 0.03 

340247.770 

CN 


3 ^ 2,7 = 7/2^ 5/2, F 

= 9/2 ^ 7/2 





340714.155 

SO 


78 ^ 67 


1.1 

0.4 

0.10 

0.04 + 0.01 

362736.048 

H 2 CO 


5o5 ^ 4o4 


1.8 

1.4 

0.23 

0.35 + 0.04 

363945.894 

H 2 CO 


524 ^ 423 


2.2 

1.5 

0.11 

0.19 + 0.04 

Table B.15. Observed line strengths 

towards SMM 2, total flux (both source and outflow components). 



Frequency 

Molecule 


Transition 


VlSR 

Av 

Tmb 

/ T^hdv 

[MHz] 





[km s"^] 

[km s"^] 

[K] 

[Kkms-‘] 

216112.582 

DCO+ 


7 = 3^2 


4.8 

0.5 

0.36 

0.21 ± 0.03 

216278.756 

C-C 3 H 2 


3o3 ^ ^21 


5.0 

0.5 

0.26 

0.14 + 0.02 

216373.320 

C 2 D 


3 ^2,7 = 7/2^5/2,F 

= 5/2 ^ 3/2 

5.4 

0.8 

0.09 

0.08 + 0.03 

216373.320 

C 2 D 


3 ^ 2,7 = 7/2^ 5/2, F 

= 7/2 ^ 5/2 





217104.919 

SiO 


7 = 5^4 


1.3 

6.0 

0.14 

1.08 + 0.08 

217238.538 

DON 


7 = 3^2 


4.3 

2.7 

0.13 

0.28 + 0.03 

217822.148 

C-C 3 H 2 


6o6 ^ 5i5 


5.0 

1.4 

0.19 

0.26 + 0.03 

217822.148 

C-C 3 H 2 


6i6 ^ 5o5 






218222.192 

H 2 CO 


3o3 ^ 2 q 2 





8.76 + 0.07 

218440.063 

CH 3 OH 


42^3i,E 


2.8 

3.9 

0.28 

1.14 + 0.05 

218475.632 

H 2 CO 


322 ^ ^21 


2.6 

3.3 

0.49 

1.77 + 0.06 

218760.066 

H 2 CO 


321 ^ ^20 


2.6 

3.2 

0.53 

1.86 + 0.05 

219560.354 

c‘«o 


7 = 2^1 


5.3 

1.3 

9.35 

14.74 + 0.07 

219908.525 

H'^CO 


3i2 ^ 2ii 





0.31+0.04 

219949.442 

SO 


65 ^ 54 





4.18 + 0.07 

337061.050 

c‘’o 


7 = 3^2 


5.3 

1.5 

1.34 

2.30 + 0.04 

338344.588 

CH 30 H 


7-1 ^ 6 - 1 , E 





1.78 + 0.07 

338408.698 

CH 30 H 


7o ^ 60 , A+ 





1.76 + 0.06 

340019.626 

CN 


3 ^ 2,7 = 5/2^ 3/2, F 

= 3/2 ^ 3/2 

4.9 

0.8 

0.09 

0.09 + 0.02 

340031.549 

CN 


3 ^2,7 = 5/2^3/2,F 

= 7/2 ^ 5/2 

5.0 

2.7 

0.39 

1.96 + 0.05 

340035.408 

CN 


3 ^ 2,7 = 5/2^ 3/2, F 

= 5/2 ^ 3/2 





340035.408 

CN 


3 ^2,7 = 5/2^3/2,F 

= 3/2 ^ 1/2 





340247.770 

CN 


3 ^ 2,7 = 7/2^ 5/2, F 

= 7/2 ^ 5/2 

4.9 

2.1 

1.12 

3.13 + 0.04 

340247.770 

CN 


3 ^2,7 = 7/2 ^5/2, F 

= 9/2 ^ 7/2 





340261.773 

CN 


3 ^ 2,7 = 7/2^ 5/2, F 

= 5/2 ^ 5/2 

5.0 

1.2 

0.09 

0.12 + 0.02 

340264.949 

CN 


3 ^2,7 = 7/2^5/2,F 

= 7/2 ^ 7/2 

4.8 

0.9 

0.09 

0 . 11 + 0.02 

340714.155 

SO 


78 ^ 67 


3.0 

3.1 

0.09 

0.25 + 0.03 

362630.303 

HNC 


7 = 4^3 


4.8 

0.9 

1.52 

1.43 + 0.05 

362736.048 

H 2 CO 


5o5 ^ 4o4 





1.63 + 0.10 
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Table B.16. Observed line strengths towards SMM 2, on-source velocity component. The integrated intensities are calculated from Gaussian fits. 


Frequency 

[MHz] 

Molecule 

Transition 

VlSR 
[km s"^] 

Av 

[km s"^] 

Tmb 

[K] 

/ T„i,dv 
[Kkms-‘] 

218222.192 

H 2 CO 

3o3 ^ 2 o 2 

5.1 

1.3 

0.87 

1.23 + 0.02 

218440.063 

CH 3 OH 

42^3i,E 

5.2 

0.8 

0.11 

0.09 + 0.02 

218475.632 

H 2 CO 

322 ^ 221 

5.1 

1.6 

0.08 

0.14 + 0.03 

218760.066 

H 2 CO 

321 ^ 220 

4.6 

1.5 

0.16 

0.26 + 0.02 

219908.525 

0 

U 

r-H 

X 

3i2 ^ 2 ii 

4.6 

2.2 

0.07 

0.16 + 0.03 

219949.442 

so 

65 ^ 54 

5.7 

0.5 

0.23 

0.13 + 0.01 


338344.588 

CH 3 OH 

7-1 ^ 6 _i,E 

6.4 

2.5 

0.04 

0 . 11 + 0.02 

338408.698 

CH 3 OH 

2o ^ 60 , A+ 

5.2 

1.6 

0.05 

0.09 + 0.02 

362736.048 

H 2 CO 

5o5 ^ 4o4 

5.3 

1.7 

0.20 

0.35 + 0.05 


Table B.17. Observed line strengths towards SMM 2, outflow velocity component. The integrated intensities are calculated from Gaussian fits. 


Frequency 

Molecule 

Transition 

vlsr 

Av 

Tmb 

/ Ttnbdv 

[MHz] 



[km s-‘] 

[km s-‘] 

[K] 

[Kkms-‘] 

217104.919 

SiO 

7 = 5^4 

1.3 

6.6 

0.32 

1.08 + 0.08 

218222.192 

H 2 CO 

3o3 ^ 2 o 2 

2.7 

2.9 

2.11 

6.55 + 0.03 

218440.063 

CH 3 OH 

42^3i,E 

2.6 

3.3 

0.29 

1.02 + 0.04 

218475.632 

H 2 CO 

322 ^ 221 

2.5 

2.8 

0.51 

1.54 + 0.03 

218760.066 

H 2 CO 

321 ^ 220 

2.4 

2.4 

0.56 

1.44 + 0.03 

219908.525 

H^^CO 

3i2 ^ 2ii 

2.0 

1.8 

0.09 

0.17 + 0.03 

219949.442 

SO 

65 ^ 54 

3.1 

3.1 

0.96 

3.20 + 0.04 

338344.588 

CH 3 OH 

7-1 ^ 6 - 1 , E 

2.6 

3.2 

0.23 

0.76 + 0.03 

338408.698 

CH 3 OH 

2o ^ 60, A+ 

2.6 

2.5 

0.32 

0.85 + 0.02 

362736.048 

H 2 CO 

5o5 ^ 4o4 

2.7 

2.9 

0.43 

1.31 + 0.06 


Table B.18. Observed line strengths towards CX042. 


Frequency 

Molecule 

Transition 

VlSR 

Av 

Tmb 

f T^dv 

[MHz] 



[km s"^] 

[km s-‘] 

[K] 

[Kkms-‘] 

218222.192 

H 2 C 0 

%3 2o2 

5.8 

1.9 

0.07 

0.14 + 0.02 

219560.354 

c‘*o 

7 = 2^1 

5.4 

1.2 

2.40 

3.34 + 0.02 

337061.050 

c‘^o 

J^3^2 

5.4 

0.8 

0.20 

0.16 + 0.02 
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Table B.19. Observed line strengths towards CrA-44. 


Frequency 

Molecule 

Transition 


VlSR 

Av 

Tmb 

f Tn,bdv 

[MHz] 




[km s"^] 

[km s"^] 

[K] 

[Kkms-*] 

216112.582 

DCO+ 

7 = 3^2 


5.7 

0.6 

0.83 

0.51+0.02 

216278.756 

C-C 3 H 2 

3o3 ^ 221 


5.7 

0.8 

0.71 

0.64 + 0.02 

216373.320 

C 2 D 

A = 3 ^ 2 ,/ = 7/2^5/2,F = 

5/2 ^ 3/2 

6.0 

1.4 

0.20 

0.28 + 0.02 

216373.320 

C 2 D 

A = 3 ^ 2 ,/ = 7/2^5/2,F = 

7/2 ^ 5/2 





216428.320 

C 2 D 

A = 3 ^ 2 ,/ = 5/2^3/2,F = 

7/2 ^ 5/2 

5.5 

1.0 

0.18 

0.19 + 0.02 

216428.320 

C 2 D 

A = 3 ^ 2 ,/ = 5/2^3/2,F = 

5/2 ^ 3/2 





217238.538 

DCN 

7 = 3^2 


5.6 

1.1 

0.19 

0.24 + 0.02 

217467.150 

'^CN 

N = 1,7 = 5/2^ 3/2, Fi = 3^2,F = 4^3 

5.3 

1.1 

0.08 

0.10 + 0.02 

217467.150 

‘^CN 

TV = 2 ^ 1,7 = 5/2 ^ 3/2,Fi = 3 - 

2, F = 3 ^ 2 





217822.148 

C-C 3 H 2 

6 o 6 ^ 5 i 5 


5.7 

1.0 

0.65 

0.73 + 0.02 

217822.148 

C-C 3 H 2 

6 i 6 ^ 5 o 5 






217940.046 

C-C 3 H 2 

5 i 4 ^ 423 


5.7 

0.8 

0.47 

0.42 + 0.02 

218160.456 

C-C 3 H 2 

524 ^ 4 i 3 


5.8 

1.2 

0.13 

0.17 + 0.02 

218222.192 

H 2 CO 

3 o 3 ^ 2 o 2 


5.7 

0.9 

1.11 

1.18 + 0.02 

219560.354 


7 = 2^1 


5.7 

0.9 

7.41 

7.29 + 0.02 

219949.442 

so 

65 ^ 54 


5.8 

1.0 

0.86 

1.06 + 0.03 

337061.050 


7 = 3^2 


5.6 

1.1 

1.30 

1.60 + 0.04 

340031.549 

CN 

A = 3^2,7 = 5/2^3/2,F = 

7/2 ^ 5/2 

5.5 

1.4 

0.33 

0.47 + 0.03 

340035.408 

CN 

A = 3 ^2,7 = 5/2^3/2,F = 

3/2 ^ 1/2 

5.6 

1.1 

0.41 

0.44 + 0.03 

340035.408 

CN 

N = 3 ^2,7 = 5/2^3/2,F = 

5/2 ^ 3/2 





340247.770 

CN 

A = 3 ^2,7 = 7/2^5/2,F = 

7/2 ^ 5/2 

5.4 

1.8 

0.82 

1.57 + 0.04 

340247.770 

CN 

N = 3 ^2,7 = 7/2^5/2,F = 

9/2 ^ 7/2 





340714.155 

SO 

78 ^ 67 


5.5 

1.5 

0.33 

0.51 + 0.03 


Table B.20. Observed line strengths towards CrA-33. 


Frequency 

[MHz] 

Molecule 

Transition 

VlSR 

[kms-‘] 

Av 

[km s-‘] 

Tmt 

[K] 

f Tynbdv 
[Kkms-‘] 

216112.582 

DCO+ 

7 = 3^2 

5.7 

0.4 

0.55 

0.26 + 0.02 

216278.756 

C-C 3 H 2 

3o3 ^ 221 

5.7 

0.1 

0.17 

0.06 + 0.02 

218222.192 

H 2 CO 

3o3 ^ 2o2 

5.7 

0.4 

0.76 

0.31 + 0.02 

218440.063 

CH 3 OH 

42^3i,E 

5.7 

0.5 

0.17 

0.10 + 0.03 

219560.354 

C^^O 

7 = 2^1 

5.7 

0.8 

4.51 

3.92 + 0.03 

219949.442 

SO 

65 ^ 54 

5.7 

0.3 

0.45 

0.15 + 0.02 

337061.050 

c‘^o 

7 = 3^2 

5.6 

0.7 

0.24 

0.18 + 0.02 


Table B.21. Observed line strengths towards VV CrA. 


Frequency 

[MHz] 

Molecule 


Transition 


VlSR 
[km s"^] 

Av 

[km s"^] 

^mb 

[K] 

f Tynbdv 
[Kkms-‘] 

216112.582 

DCO+ 


7 = 3^2 


4.9 

0.7 

0.15 

0.12 + 0.01 

216278.756 

C-C 3 H 2 


3o3 ^ 221 


5.0 

0.7 

0.13 

0.09 + 0.01 

217238.538 

DCN 


7 = 3^2 


4.8 

0.8 

0.05 

0.04 + 0.01 

217822.148 

C-C 3 H 2 


6o6 ^ 5i5 


5.0 

0.8 

0.13 

0.12 + 0.01 

217822.148 

C-C 3 H 2 


6i6 ^ 5o5 






217940.046 

C-C 3 H 2 


5i4 ^ 423 


5.1 

0.8 

0.08 

0.07 + 0.01 

218222.192 

H 2 CO 


3o3 ^ 2o2 


4.9 

0.9 

0.78 

0.79 + 0.02 

218475.632 

H 2 CO 


322 ^ 221 


4.8 

1.3 

0.06 

0.09 + 0.02 

218760.066 

H 2 CO 


321 ^ 220 


5.0 

1.4 

0.07 

0.12 + 0.02 

219560.354 

C^^O 


7 = 2^1 


4.9 

0.9 

5.48 

6.12 + 0.02 

219949.442 

so 


65 ^ 54 


5.0 

1.3 

0.84 

1.39 + 0.02 

337061.050 

c^^o 


7 = 3^2 


4.9 

1.0 

0.76 

0.86 + 0.02 

340031.549 

CN 

N = 

3 ^2,7 = 5/2^3/2,F 

= 7/2 ^ 5/2 

4.8 

1.4 

0.20 

0.29 + 0.02 

340035.408 

CN 


3 ^ 2,7 = 5/2^ 3/2, F 

= 3/2 ^ 1/2 

4.8 

1.2 

0.19 

0.25 + 0.02 

340035.408 

CN 


3 ^2,7 = 5/2^3/2,F 

= 5/2 ^ 3/2 





340247.770 

CN 


3 ^ 2,7 = 7/2^ 5/2, F 

= 7/2 ^ 5/2 

4.7 

1.5 

0.57 

0.95 + 0.02 

340247.770 

CN 


3 ^2,7 = 7/2 ^5/2, F 

= 9/2 ^ 7/2 





340714.155 

SO 


78 ^ 67 


4.8 

2.0 

0.31 

0.69 + 0.03 
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Table B.22. Observed line strengths towards CrA-37. 


Frequency 

Molecule Transition 

Vlsr 

Av 

^mb 

f 7’mbrfv 

[MHz] 


[km s"^] 

[km s"^] 

[K] 

[Kkms-*] 

219560.354 

T 

(N 

II 

o 

00 

u 

5.7 

0.6 

1.12 

0.72 ± 0.03 
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Appendix C: Complete IRS7B spectrum at 217-245.5 GHz 

All spectra in this and the next section are cut off at 0.5 K to also show the weaker spectral lines. 
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Fig. C.2. See Fig. CL 
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Fig. C.3. See Fig. CL 
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Fig. CA See Fig. CL 
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Fig. C.7. See Fig. CL 
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Appendix D: Complete spectra of the H 2 CO, CH 3 OH, and C-C 3 H 2 survey 



Fig. D.l. 0.9 mm spectrum of IRS7B, smoothed by a factor 8, corresponding to a channel width of 0.6 km 
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Rest Frequency [GHz] 

CrA-3 



Rest Frequency [GHz] 



Rest Frequency [GHz] 

Fig. DA 1.4 mm spectrum of CrA-3, smoothed by a factor 8, corresponding to a channel width of 0.8 km s ^ 
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CrA-3 
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0.3 - 


0.2 
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Fig.D.5. 0.9 mm spectrum of CrA-3, smoothed by a factor 8, corresponding to a channel width of 0.6 km s ^ 
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Rest Frequency [GHz] 



Rest Frequency [GHz] 

Fig.D.6. 1.4 mm spectrum of CrA-5, smoothed by a factor 8, corresponding to a channel width of 0.8 km s 
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LS-RCrAl 
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Fig.D.9. 1.4 mm spectrum of LS-RCrAl, smoothed by a factor 8, corresponding to a channel width of 0.8 km s ^ 
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Fig.D.ll. 1.4 mm spectrum of Haas4, smoothed by a factor 8, corresponding to a channel width of 0.8 km s ^ 
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Fig.D.16. 0.9 mm spectrum of IRS5A, smoothed by a factor 8, corresponding to a channel width of 0.6 km s ^ 
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Fig.D.18. 1.4 mm spectrum of IRS5N, smoothed by a factor 8, corresponding to a channel width of 0.8 km s ^ 
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Fig.D.21. 1.4 mm spectrum of IRSl, smoothed by a factor 8, corresponding to a channel width of 0.8 km s ^ 
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IRS7A 
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Fig. D.25. 0.8 mm spectrum of IRS7A, smoothed by a factor 8, corresponding to a channel width of 0.5 km s ^ Strong atmospheric lines are 

visible at 364.1 GHz, 364.4 GHz, and 365.3 GHz. 
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CrA-24 
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mm spectrum of CrA-24, smoothed by a factor 8, corresponding to a channel width of 0.5 km s 
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SMM 2 
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Fig.D.29. 1.4 mm spectrum of SMM 2, smoothed by a factor 8, corresponding to a channel width of 0.8 km s ^ 
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SMM 2 
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Fig. D.30. 
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Fig. D.32. 
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CX042 
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1.4 mm spectrum of CX042, smoothed by a factor 8, corresponding to a channel width of 0.8 km s“^ 
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Fig. D.33. 
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W CrA 
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Appendix E: Unidentified spectrai iines 

We here present our attempts to identify all U lines detected to¬ 
wards IRS 7 B: The faint line at 229.7760 GHz (assuming vlsr = 
5.7 km s-i) could be the CH3CHO line at 229.7750 GHz, but that 
would give this line a vlsr = 4.5 km s“\ which is 1.2 km s“^ 
lower than the median LSR velocity and 0.5 km s“^ lower than 
any other line, including the more certain CH3CHO lines (also 
if taking the JPL uncertainty of 50 kHz on the measured labo¬ 
ratory line frequency into account). The line at 225.160 GHz 
coincides with a CH3SH line, but this is not the strongest line 
expected from this species in the studied frequency range. The 
237.998 GHz line is consistent with a c-H^^CCCH line, but if 
this is a true detection, the c-H^^CCCH abundance would be 
of the same order as the C-C3H2 abundance, which is very un¬ 
likely. The remaining U lines do not align with any spectral 
lines found in Splatalogue of species expected to be present at 
this rms level. Two of the detected lines agree with known U 
lines at 223.756 GHz and 233.456 GHz. As a result of incom¬ 
plete image-band rejection, a few so-called ghost lines appear 
at the image frequencies of the strongest spectral lines (CO iso- 
topologue lines and strong H2CO, CH3OH, and CN lines). We 
did not tabulate these ghost lines, but they are labelled in the 
spectra in Appendix C. For all identified and unidentified lines, 
the image frequencies were investigated to exclude the possibil¬ 
ity that any reported line would be such a ghost line. 

Appendix F: Non-LTE models of H2CO emission 

Mangum & Wootten ( 1993 ) showed that ratios of H2CO tran¬ 
sitions involving the same /^-level but from different ^-ladders 
(such as the 3o3 ^ 2o2/322 ^ 221 and 5o5 ^ 4o4/524 ^ 423 
ratios) are excellent tracers of the kinetic temperature of the gas 
because they only operate through collisional excitation. How¬ 
ever, at densities ^ < 10^ cm“^, the different J^-lewels are not 
fully thermalised, which means that ratios of transitions involv¬ 
ing different /u-levels (such as 3o3 ^ 2o2/5o5 ^ 4o4) are sensi¬ 
tive not only to the temperature, but also to the molecular density 
^(H2). 

Mangum & Wootten ( 1993 ) used LVG models to derive a 
method for extracting temperature, density, and column den¬ 
sity from H2CO line observations. Jansen ( 1995 ) performed 
RADEX modelling of H2CO line ratios to show which H2CO 
line ratios can be used to trace the kinetic temperature at rather 
low column densities, investigating certain line ratios at a p- 
H2CO column density of A = 10 ^^ cm“^. In Fig. F.l we 
show such plots at a P-H2CO column density N = 10 ^^ cm“^, 
which agrees better with the properties of the sources in this 
study. The 3o3 ^ 202/5o5 ^ 4o4 ratio clearly is a particu¬ 
larly poor temperature probe, while the 3o3 ^ 2o2/322 ^ 221, 
3 o3 ^ 2o2/32i ^ 220, and 5o5 ^ 4o4 / 523 ^ 422 ratios probe the 
temperature well at T < 50 K and ^ > 10 ^ cm“^, also at this rela¬ 
tively high column density. We used the P-H2CO molecular data 
file from the LAMDA database (Schoier et al. 2005 ), which uses 
the collisional rates for P-H2CO-H2 from Wiesenfeld & Faure 
( 2013 ). If we instead use the older Green ( 1991 ) p-H2CO-He 
collisional rates corrected to collisions with H2 by correctional 
factors for pressure broadening and relative collision velocities 
the estimated ^(H2) values become ~ 50 % higher. 

To acquire the H2CO temperature T, column density A, and 
the H2 density n towards the sources where we have measure¬ 
ments of both Ju = 2 and Ju = ^ transitions, we used the follow¬ 
ing recipe developed from the method of Mangum & Wootten 
( 1993 ): 


1 . Measure the rotational temperatures of the Ju = 2 and Ju = ^ 
transitions separately and assume that the kinetic tempera¬ 
ture is the weighted average of these two temperatures. If 
only the = 3 transitions can be measured, we use this 
as the kinetic temperature (and conversely, we use only the 
/u = 5 temperature in IRS 7 A, since the Ju = 2 transitions 
were not covered). 

2 . Use the non-LTE radiative transfer code RADEX (van der 
Tak et al. 2007 ) to calculate model line strengths of the 3o3 ^ 
2 o 2 and 5o5 ^ 4o4 lines for the kinetic temperature and a 
grid of n and A (where A is the total H2CO column density 
assuming an ortho-to-para ratio of 1.6). 

3 . On this (n,N) model grid, plot lines corresponding to the 
measured 3o3 ^ 2o2/5o5 ^ 4o4 line ratio and the measured 
3 o3 ^ 2o 2 line strength. The values of n and A are deter¬ 
mined from the intersection of these lines (see example in 
Fig. F. 2 ). 
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Fig. E2. RADEX fit for H 2 CO in IRS7B, given 7 = 40 K (from the 
rotational diagram fit of the full IRS7B survey). The solid line shows 
the measured 3o3 ^ 2 o 2 line intensity, and the dashed line is the mea¬ 
sured 3 o3 ^ 2o2/5o5 ^ 4o 4 line ratio. 3cr errors are shown as grey 
lines, but can barely be distinguished from the measured values. When 
we estimate uncertainties, the errors on the fitted temperatures are also 
taken into account. 
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